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SUMMARY 


This  program  is  concerned  with  the  development  of  an  interesting  class  of  artificial 
dielectrics  and  a  deeper  understanding  of  its  electromagnetic  properties  at  high  frequencies 
(centimeter  and  millimeter  wavelengths).  This  Final  Technical  Report-Part  A  deals  with  essentially 
non-magnetic  composite  dielectrics.  Part  B,  which  contains  classified  material,  treats  magneto¬ 
dielectric  composites. 

A  model  dielectric  consisting  of  a  randomly  packed,  metal-insulator  composite  is  an 
example  of  an  inhomogeneous  material  whose  properties  have  provided  problems  and  challenges 
for  theorists  and  experimentalists  for  over  150  years.  We  have  studied  the  permeability  and 
permittivity  at  high  frequencies  on  correlated  artificial  dielectrics  with  moderately  sized  metal 
particles  (-1  |im  to  -37  (im)  randomly  loaded  into  both  polymeric  and  inorganic  glassy  binders  up 
to  40  vol%.  The  metal  alloy  powder,  Ni92Cr8,  is  prepared  by  gas-water  atomization,  which 
provides  an  oxide  coating  for  interparticle  isolation  (■  the  correlation).  The  powder  is 
ferromagnetic  below  about  160K,  which  enables  a  low-temperature  magnetic  measurement  of 
volume  loading,  p,  while  leaving  a  “non-magnetic"  composite  at  room  temperature.  Non¬ 
percolating  composites  with  p  ranging  from  0.1  to  0.4,  for  separate  size  fractions  (in  |im)  of -10, 
+10-20,  +20-37,  were  examined  from  about  1  GHz  to  20  GHz  and  for  selected  samples  at  -35 
GHz  and -90  GHz. 

The  complex  permeability  results  are  in  very  good  agreement  with  induced  magnetic  dipole 
(eddy  current)  calculations  as  a  function  of  frequency  and  particle  diameter,  Permittivity  results 
(real  pan)  are  independent  of  particle  diameter  and  essentially  constant  with  frequency.  Their 
dependence  on  loading  is  compared  to  the  classical  Maxwell/Qausius-Mossotti  dipolar  calculation, 
to  the  ordered  array  models  with  higher-order  multipole  interactions  pioneered  by  Rayleigh  and 
completed  recently  by  Doyle  and  McPhedran  et  al.,  as  well  as  to  a  pair-interaction  model  developed 
during  this  program.  Comparison  is  also  made  to  an  empirical  curve  obtained  from  remarkably 
ideal  experiments  of  Guillien  and  of  Turner.  Our  permittivity  data  on  another  alloy  powder 
system,  from  work  outside  this  contract,  are  in  excellent  agreement  with  the  empirical  curve.  For 
the  Ni-Cr  alloy  powder  composites,  there  is  a  residual  upward  deviation  of  ef  vs.  p,  which  is 
attributed  to  (detectable)  non-spherical  deviations  in  particle  morphology. 

This  research  demonstrates  that  the  microwave  electromagnetic  behavior  of  non¬ 
percolating,  artificial-dielectric  composites  is  well  understood.  Random  arrays  have  their  own 
unique  properties,  as  distinct  from  ordered  arrays.  The  present  study  provides  a  firm  foundation 
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I.  INTRODUCTION 


This  research  is  concerned  with  the  development  of  a  novel  class  of  artificial  dielectrics  and 
a  deeper  understanding  of  its  properties  at  millimeter  wavelengths.  Our  original  point  of  departure 
was  a  Radar  Absorbing  Materials  Physics  Model  developed  during  a  classified  program  [1].  That 
study  involved  experimental  measurements  in  the  microwave  range  of  frequencies  (centimeter 
wavelength)  over  a  wide  range  of  parameter  space  and  analysis  of  the  broad  data  base  obtained 
therefrom.  Many  of  the  model  features  and  predictions  are,  however,  relevant  to  the  millimeter 
wavelength  region.  Some  of  the  interesting  implications  of  the  prior  work  concern  the  permittivity 
behavior  of  the  composite  system,  one  which  may  be  regarded  as  a  new  class  of  oxide-passivated 
metal-insulator  dielectrics. 

A  somewhat  optimistic  overview  of  the  research  goals  for  the  study  was  stated  in  its  early 
stages  in  question  form  as  follows:  How  valid  is  the  RAM  physics  model  in  higher  frequency 
regimes?  Why  does  the  dielectric  behavior  appear  to  follow  the  simple  Maxwell  approximation, 
given  the  pitfalls  that  often  beset  metal-insulator  composites?  Has  the  artificial  (magneto-) 
dielectric  system  been  optimized  in  the  centimeter  wavelength  regime?  Can  composites  be 
fabricated  which  work  well  in  the  millimeter  wavelength  regime?  How  do  they  compare  with  other 
materials  approaches? 

Alternatively,  a  more  formal  view  of  the  program  can  be  obtained  from  the  Statement  of 
Work  given  in  the  proposal  and  reproduced  below  as  Table  1.  Concurrent  with  this,  we  take  note 
of  the  division  of  the  work  and  of  die  reports  into  two  subclasses,  one  of  essentially  non-magnetic 
composite  dielectrics  and  one  of  composite  magneto-dielectrics.  Some  aspects  of  the  work  on  the 
latter  are  classified.  Therefore  this  Final  Report  is  prepared  in  two  sections,  Parts  A  and  B,  with 
the  magneto-dielectrics  work  appearing  in  the  classified  section,  Part  B.  The  Statement  of  Work 
must  be  considered  from  that  divided  point  of  view. 

The  notation  in  the  table  following  Task  4  signals  a  modification  approved  by  the  ONR 
Scientific  Officer  (W.R.  Schmidt,  CDR,  USN)  during  an  oral  review/site  visit  on  12  June  1984.  It 
resulted  from  the  early  realization  that  particle  size  effects  and  induced  eddy  currents  would  be 
rather  important  and  that  these  had  not  been  planned  for  study  in  the  original  program  plans.  An 
expansion  of  scope  in  Tasks  3  and  6  to  these  ends-i.e.,  theoretically  and  experimentally-formed  a 
balance  to  the  suppression  of  Task  4. 
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Tabic  1 


STATEMENT  OF  WORK  BY  TASK 
TASK 

1.  Investigate  effect  of  binder  permittivity  in  artificial  dielectrics 

2.  Analyze  effects  of  microstructure  control  in  passivated  metal-insulator  composites 

3.  Develop  alternative  theoretical  approaches  to  composite  dielectrics 

4.  Search  for  percolation  transition:  exploitation  and  control 

5.  Explore  alternative  passivation  coatings/investigate  high-frequency  properties 

6.  Quantify  skin-depth  effects  in  passivated  artificial  dielectrics 

7.  Investigate  millimeter-wavelength  properties  of  artificial  dielectric  composites 

8.  Investigate  or  simulate  temperature  dependence  of  high-frequency  properties 


SUBCLASS 

ASSIGNMENT 

A 

A 

A 

removed 
see  note  in  text 

B 

A  B 
A  B 

B 
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The  program  plan  originally  provided  for  a  three-year  effort,  starting  officially  on  1  Aug. 
1983.  The  studies  under  Part  A  were  completed  by  30  June  1987  while  those  of  Part  B  continued 
into  the  latter  part  of  1988.  These  were  “no-cost”  extensions.  Reasons  for  the  delays  in  research 
and  in  reporting  were  several:  research  leaves  by  the  principal  investigator,  measuring  equipment 
problems,  and  considerable  difficulties  in  the  preparation  and/or  acquisition  of  desired  materials. 
We  regret  these  substantial  delays,  However,  we  believe  that  we  have  completed  the  tasks  of  the 
work  statement  and  contributed  to  understanding  and  progress  in  a  number  of  the  designated  areas 
of  research.  We  trust  the  reader  will  concur  after  examination  of  these  final  reports. 

The  principal  investigator  under  this  contract  has  benefitted  from  major  and  minor 
contributions  from  a  number  of  collaborators  and  consultants.  It  is  a  pleasure  to  acknowledge  their 
help,  as  listed  here  with  their  areas  of  contributions; 

L. M.  Levinson,  for  original  program  planning  and  helpful  counsel 

W.O.  Moffatt  for  alloy  preparation 

S.A.  Miller  with  D.N.  Wernple  for  alloy  powder  preparation  by  gas-atomization  and 
exploration  of  alternative  coatings 

A.E.  Berkowitz  and  J.L.  Walter  with  J.J.  Rogers  and  W.E.  Rollins  for  experimental 
research  on  the  spark  erosion  method  of  alloy  powder  preparation 

R. W.  Harrison,  D.B.  Lindeman,  and  J.  Murphy  for  consultation  on  and  preparation  of 

polymeric  composites  (OE  Aircraft  Engine,  Cincinnati); 

M. P.  Borom  with  L.E.  Szala  for  searching  for  alternative  inorganic  binders  and  preparing 
composites  thereof 

S. P.  Yu  and  J.P.  Quine  for  consultation  on  microwave  measurement  techniques 

W.B.  Westphal,  Prof.  F.R.  Morgenthaler,  and  E.  Maxwell  (MIT,  during  leave  9/83-1/84) 
for  discussions  on  high  frequency  measurements  and  on  artificial  dielectrics 

J.Q.  Hanson  for  microwave  measurements,  Ka  band,  (GE-Re-Entry  Systems,  Philadelphia) 

D.  W.  Forester  and  F.J.  Rachford  for  consultation  and  measurements  both  at  centimeter  and 
millimeter  wavelengths  (Naval  Research  Laboratory) 

D.R.  Kelland  for  consultation  on  fine  particle  separation  (National  Magnet  Lab.,  MIT) 

H.J.  Patchen  for  a  wide  range  of  important  measurement  responsibility  for  the  principal 
investigator  including  thermomagnetic,  microwave  and  Mossbauer  spectroscopy 

Prof.  W.T.  Doyle  for  major  contributions  in  the  modeling  of  random  artificial  dielectrics 
(Dartmouth) 

Prof.  W.P.  Wolf  for  a  wide  range  of  consultations  on  physics  and  measurements  (Yale). 
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II.  DEVELOPMENT  OF  AN  ARTIFICIAL  DIELECTRIC 


With  the  initial  emphasis  of  the  program  on  the  dielectric  aspects,  an  important  phase  has 
been  the  development, of  a  suitable  laboratory  artificial  dielectric.  The  model  of  a  dielectric  which 
consists  of  a  series  of  metallic  globules  separated  (Tom  each  other  by  insulating  material  goes  back 
at  least  to  Faraday  in  1837.  An  excellent  history  of  inhomogeneous  media  concepts  has  been 
recorded  by  Landauer  [2]. 


II.I  Criteria 

Our  criteria  were  that  the  metal  (alloy)  powder  be  non-magnetic  at  room  temperature  but 
ferromagnetic  at  low  temperature  to  facilitate  determination  of  the  volume  loading,  p.  The  latter 
quantity  is  key  to  quantifying  the  dielectric  behavior.  We  also  required  that  the  powder  be  capable 
of  developing  a  thin  insulating  oxide  layer  on  each  particle.  We  contrast  this  with  some  other  work 
in  which  researchers  have  encouraged  formation  of  a  thin  oxide  coating  in  powder  preparation  to 
discourage  initial  agglomeration  but  kept  it  thin  enough  to  allow  metal-to-metal  contract  between 
the  particles  under  sufficient  pressure  [3]. 

The  choice  of  candidate  alloys  thus  involved  magnetic  and  metallurgical  consideration. 
Added  to  these  was  the  challenge  of  fine  metal  powder  preparation.  For  this  task  we  made  use  of 
our  gas-atomization  facilities.  In  this  process  [4]  molten  alloy  is  forced  through  a  specially 
designed  nozzle,  and  the  liquid  stream  is  disrupted  by  high-velocity  gas  blasts  to  create  fine 
droplets.  These  solidify  as  they  descend  in  a  large  volume  settling  chamber.  In  one  of  our 
facilities  liquid-substrate  quenching  is  a  special  feature;  i.e.,  a  blast  of  water  spray  follows  the  gas 
jets  (see  Figure  1). 


II.2  Consideration  of  Candidate  Metallic  Systems 

Desirable  alloy  systems  were  therefore  such  that  one  of  the  elemental  components  would 
tend  to  form  a  stable  and  at  least  partially  protective  oxide.  Some  usual  examples  of  this  category 
are  Al,  Cr.  and  Y.  Folding  these  in  with  the  magnetic  considerations  yielded  several  suggestions, 
e.g.,  Ni  with  Cr,  Od  with  Y,  and  (Cu-Ni)  with  Al.  The  first  two  systems  are  reasonably  well 
known  magnetically  in  handbooks  [5]  or  in  more  recent  research  [6],  while  the  third  has  had  little 
exploration  [7].  The  Gd-Y  combination  was  quickly  rejected  (by  our  metallurgical  colleagues)  as 
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explosively  undesirable.  Although  the  solubility  of  A1  in  the  Ni-Cu  host  is  modest,  at  about  10 
at%  (-5  wt%),  there  appeared  to  be  a  fair  prospect  of  a  Curie  point  below  room  temperature.  We 
chose  to  prepare  an  exploratory  ingot  with  the  composition  (by  weight)  Ni90Cu5A15.  Finally  for 
Ni-Cr,  the  wide  range, of  solubility  and  well-known  magnetic  properties  enabled  us  to  specify 
“with  confidence"  an  ingot  composition  of  Ni92Cr8  with  every  expectation  of  a  Curie  point  near 
240K,  and  a  low  temperature  saturation  magnetization  of  about  <7M  (OK)  ■  22  emu/g. 


11,3  Magnetic  Characterization  of  Candidate  Metallic  Systems 

Thermomagnetic  scans,  magnetization  vs.  temperature  at  50  Oe,  along  with  low- 
temperature  (6K,  20K)  saturation  magnetization  data,  were  measured.  For  this  purpose  we  have  a 
vibrating  sample  magnetometer  system  with  capabilities  from  4K  to  1370K  (~1000*C)  and  fields 
up  to  25  kOe.  The  scans  for  Ni-Cr  are  shown  in  Figure  2.  As  for  the  Ni-Cu-Al  ternary  (not 
shown),  its  Tc  centered  at  256K,  but  there  was  a  risk  of  some  residual  undesirable  ferromagnetism 
at  room  temperature,  The  opposite  problem  presented  itself  with  the  ingot  Ni92Cr8,  Both  the 
Curie  point  and  the  magnetization  fell  short  of  confident  expectations  by  a  factor  of  2!  In  searching 
for  the  "fault"  here,  we  verified  the  chemical  composition  of  the  ingot,  and  also  prepared  a  second 
ingot  with  composition  Ni94Cr6.  Its  Tc,  measured  as  237K,  again  compares  poorly  with  the 
literature  value  of  335K1  At  this  point  we  were  facing  an  unexpected,  annoying,  and  delaying 
problem.  A  resolution  of  the  dilemma  was  clearly  warranted  because  this  Ni-Cr  system  remained 
the  most  attracdve  of  the  ones  considered  as  candidates. 


II.3.1  Exploration  of  Magnetic  Characterization  Problems 

The  solution  of  the  problem  and  the  "validation"  of  our  own  work  came  from  a  detailed 
look  at  the  published  literature  both  backward  from  the  handbook  data  and  forward  to  the  present 
time  via  a  library  bibliographic  search.  The  magnetic  measurements  cited  in  the  standard  references 
(5)  are  based  on  experimental  work  firom  the  Strasbourg  laboratory  directed  by  P.  Weiss,  a  world- 
famous  organization  that  set  most  of  the  standards  for  current  magnetics  research.  In  studying  the 
primary  sources,  the  research  of  Marian  [8]  and  Sadron  [9],  we  found  that  they  used  identical 
alloy  samples  which  had  been  prepared  earlier  by  S  afrdnek  [  10],  also  at  Strasbourg.  In  the 
preparation,  he  had  been  concerned  with  the  purity  of  the  available  starting  materials,  especially  the 
chromium.  Indeed,  he  ultimately  prepared  his  own  Ci  by  an  electrolytic  process,  to  combine  with 
high-purity  Mond  nickel.  Satisfactory  metallurgical  micrographs  were  shown,  indicating  single- 
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phase  alloys.  Inasmuch  as  our  Curie  points  were  markedly  lower,  a  suspicious  investigator  may 
still  want  to  raise  the  purity  question.  How  much  residual  Fe  or  Co  in  either  component  would  be 
needed  to  shift  the  results  into  "accord"?  More  tolerance  is  in  order  with  regard  to  the  magnetic 
saturation  because  the,  French  workers  had  to  extrapolate  downward  from  1 10K  toward  absolute 
zero.  However,  their  results  are  40%  to  100%  greater  than  ours! 

'  As  one  looks  forward  from  the  standard  references  into  the  most  recent  two  decades,  a 
different  picture  emerges.  A  group  at  Leiden  [1 1]  in  1962  noticed  some  differences  in  saturation 
magnetization  for  alloys  of  S.6  at%  and  9. 1  at%,  and  in  an  estimate  of  Tc  for  the  former  alloy,  but 
politely  called  these  "small  discrepancies."  Next,  a  rather  thorough  investigation  in  1972  at 
Strasbourg  by  Besnus  et  al.  [12]  (on  newly  prepared  samples  and  including  low-temperatures) 
presented  results  which  differed  quite  a  bit  from  the  early  data  though  not  quite  as  much  as  the 
Leiden  results.  They  diplomatically  omitted  to  mention  that  they  differed  from  their  colleagues  of 
an  earlier  generation  1  The  emerging  picture  is  displayed  in  Figure  3  for  Curie  points  and  in  Figure 
4  for  saturation  magnetization.  (Compositions  are  in  atom  percent,  as  used  by  these  research 
workers.  In  the  region  of  interest  to  us,  i.e.,  6-8  wt%,  the  amount  in  at%  is  about  one  integer 
higher).  More  recently  a  study  from  Japan  [13]  gave  results  in  close  agreement  with  Besnus  et  al., 
and  did  make  a  significant  graphical  display  of  the  differences  with  the  early  work.  Lastly, 
Simpson  and  Smith  in  1982  reported  especially  on  compositions  with  quite  low  Curie  points  [14], 
in  close  agreement  with  our  measurements. 

While  examination  of  the  comparative  results  in  Figures  3  and  4  does  not  fully  settle  the 
question  of  which  set  is  correct,  the  existence  of  a  major  discrepancy  with  the  early  work  is 
established.  The  present  work  on  the  several  ingou  sets  the  base  line  for  subsequent  magnetic 
evaluation  of  powder  properties  and  ultimately,  for  our  original  goal,  of  volume  loadings  in 
artificial  dielectric  composites. 


II.4  Preparation  and  Magnetic  Characterization  of  Metal  Powder 

At  this  stage  we  were  ready  for  powder  preparation.  We  anticipated  (correctly)  that  there 
would  be  some  composition  shift  in  this  step  toward  a  higher  Curie  point,  and  thus  we  selected  the 
ingot  Ni92Cr8  (about  9  at%  Cr)  in  order  to  keep  the  alloy  non-magnetic  at  room  temperature. 

Two  atomization  trials  of  the  N192-Cr8  alloy  were  performed.  In  both  cases  the  melt  was 
atomized  while  at  a  temperature  of  1870K  by  argon  gas  at  4.28  MPa.  The  high  melt  temperature 


8 


{*10 


ATOM  %  Cr  IN  Ni 


Figure  3.  Ferromagnetic  Curie  points  ai  a  function  of  composition  In  NI*Cr  solid  solution 
alloys  Including  work  of  other  Investigators  as  well  as  present  results 
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4.  Saturation  magnetisation  par  avaragt  atom  vs.  composition  for  NI*Cr  alloys 
Including  available  literature 


Figure  5.  Particle  slie  distributions  from  two  gas  atomisation  runs  on  Nl92Crg  alloy 


was  employed  to  reduce  the  melt  surface  tension  and  melt  viscosity  while  the  high  gas  pressure 
resulted  in  a  high  gas  jet  velocity.  (Mach  Number  -3,2).  Previous  analysis  has  shown  that  fine 
powder  sizes  are  favored  by  increased  gas  velocity  and  decreased  melt  surface  tension  and  melt 
viscosity. 

The  two  trials  were  performed  using  two  different  versions  of  the  GE-CRD-developed 
close-coupled  gas  nozzle.  The  second  run  used  an  improved  form  of  the  close-coupled  concept, 
designed  to  provide  increased  yields  of  fine  powder. 

The  results  of  the  two  runs  are  shown  in  Figure  5.  The  particle  size  distributions  between 
37  fx  m  and  250  /x  m  were  determined  by  sieving.  The  panicle  size  distributions  below  37  fx  m 
were  obtained  by  extrapolation,  assuming  the  curves  to  be  log-normal.  (This  assumption  has  been 
verified  in  other  work). 

The  thermotnagnetic  scan  of  this  powder  (sieved  to  <37  jx  m)  is  shown  in  Figure  2  as 
powder  RS67.  An  average  Curie  point  of  165K  results,  some  68K  upward  from  the  ingot,  It  is 
broader,  suggesting  a  range  of  Tc  values.  The  shift  in  Tc  is  presumed  to  arise  principally  from  the 
effects  of  the  blast  of  water  spray. 

Referring  to  Figure  3,  this  Tc  translates  in  composition  to  about  8  at%  Cr  or  7. 15  ±  0.2 
wt%  Cr.  Using  this  composition,  there  is  satisfactory  accord  of  the  saturation  magnetization 
"line”in  Figure  4  with  the  “present  work”  data. 

For  the  practical  record,  the  magnetization  values  at  low  temperature  (6K  or  20K), 
extrapolated  to  Infinite  field,  are  21.2  emu/g  and  12.0  emu/g  for  the  ingots  Ni94Cr6  and  Ni92Cr8 
respectively.  For  the  RS67  powder  we  find  16.2  emu/g,  and  for  the  RS98  powder  prepared  later 
the  value  is  14.6  emu/g. 

II.5  Evaluation  of  Surface  Condition  of  Metal  Powder 

The  shift  in  composition  during  powder  preparation  corresponds  to  a  net  depletion  of  Cr  in 
the  alloy  particle.  It  is  tempting  to  attribute  this  to  an  oxide  coating  of  Cr203.  It  is  much  more 
difficult  to  examine  the  surface  coating  of  spherical  particles  by  microprobe  techniques  than  the 
surface  scales  on  bulk  flat  samples.  Several,  moderately  successful  attempts  to  do  this  on  these 
powder  spheres  are  reported  herein.  Conventional  oxidation  studies  on  the  Ni-Cr  system,  albeit 


11 


under  isothermal  conditions,  indicate  a  moderately  complex  scaling  with  NIC,  Cr203,  Ni-Cr  spinel 
oxide  and  other  components  [15]. 


If  we  naively  (and  tentatively)  hold  to  a  view  that  the  coating  may  be  Cr203,  we  can 
estimate  its  thickness  from  a  mass  balance,  Using  a  mean  diameter  of  about  23  /j.  m  we  calculate  a 
shell  coating  whose  thickness  is  about  0.1  m,  which  is  not  unreasonable.  A  sketch  of  this  rough 
calculation  is  as  follows: 

[(w  /  o  Cr  in  ingot)  x  vol.  sphere  x  p  alloy]  ■  [(w  /  o  Cr  in  powder)  x  vol.  sphere  x  p  alloy] 

+[(w  /  o  Cr  in  oxide)  x  (vol.  oxide  shell)  x  p  oxide] 

vol.  sphere  •  4rrr3  /  3;  vol.  shell  ■  47rr*  A  r 

A  r  hi  f  / 0  Cf  fagot  •  w  /  o  Cr  in  powder)  x  (diam.)  x  p  alloy  | 

[  (w/oCrinoxide)x3x2xpoxide  J 


Let  us  examine  our  first  attempt  with  microprobe  techniques  on  the  Ni-Cr  powder.  As- 
prepared  powder  was  pressed  onto  and  into  an  aluminum  substrate  block  for  examination  by  Auger 
electron  spectroscopy.  The  count  rates  for  the  various  elements,  Ni,  Cr  and  O,  were  monitored  as 
a  funcdon  of  sputter  time.  The  latter  is  proportional  to  depth  into  the  particle  being  examined,  as 
material  is  sputtered  away.  The  probe  spot  used  was  smaller  than  the  particle  diameter  (from  the 
powder  size  fraction  +20,  -  27  p  m).  Count  rates  are  expressed  in  atomic  percent,  after  appropriate 
scaling  in  the  instrument  (Physical  Electronics  Model  600  Auger  Microprobe).  The  "depth’’ 
probing  may  not  be  completely  simple  because  the  pressure  used  in  preparing  the  specimen  may 
have  deformed  the  spherical  alloy  particle.  As  an  alternative,  we  hoped  to  obtain  a  profile  on  a 
polished  particle  as  the  probe  traversed  the  edge  coating. 

The  depth  profile  by  Auger  spectroscopy  is  shown  in  Figure  6.  While  not  the  simple  naive 
picture  we  hypothesized,  it  does  have  features  one  might  expect.  At  the  surface  (r  ■  0,  sputter 
time)  there  is  a  major  build-up  of  oxygen,  from  a  residual  low  level  representing  a  sort  of 
background;  and  a  modest,  ~50%,  build-up  of  chromium  from  about  12  at%  to  about  18  at%. 
Concurrently,  there  is  a  major  decrease  of  nickel,  albeit  not  to  the  near-zero  level  hypothesized. 

We  can  attempt  to  rationalize  these  data  in  terms  of  NiO  and  Cr2C>3,  there  being  no  evidence  for  Ni 
metal  in  the  magnetic  data.  (Oxidation  studies  of  this  system  [15]  always  assume  that  Cr  goes  to 


12 


SPUTTER  TIME  (min.) 


Figure  (.  “Depth”  profile  by  Auger  electron  spectroscopy  of  Ni-Cr  powder  compact: 
composition  vs.  sputter  time  (proportional  to  depth) 
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Cr203,  or  more  complex  spinels  or  solid  solution).  The  "surface  composition,"  Nio.32  Cro.is 
Oo,50»  is  roughly  equivalent  to  3NiO  •  O2O3,  which  corresponds  to  Nio.27Cro.i8  Oo,54- 

If  we  oversimplify  the  surface  as  a  uniform  shell  (again),  we  can  invoke  a  mass  or  atom 
balance  to  determine  how  much  material  is  in  the  shell.  Using  the  chemically  and  magnetically 
determined  alloy  compositions,  we  write: 

N10.91  Cro.09  +  6rO->  y  Nio,92  Cro.08  +  2  (3NiO-l  Cr203). 

Balance  equations  for  Ni  and  Cr  yield  two  simultaneous  equations  in  y  and  z,  which  are  readily 
solved  to  yield  y  «  0.969,  r  *  0.00625.  Approximately  2%  of  the  original  nickel  goes  into  the 
oxide  shell,  while  about  14%  of  the  original  chromium  makes  the  transition.  This  does  confirm  the 
general  stronger  tendency  for  Cr  to  oxidize.  Expressing  the  chemical  balance  in  terms  of  the  mass 
of  one  average  (23  fi  m)  particle  and  the  number  of  moles  of  alloy  in  that  panicle,  we  can  obtain  the 
masses  of  NiO  and  Cr203.  Converting  to  volumes  (via  their  separate  specific  gravity  values),  and 
equating  their  total  volume  to  a  "shell"  volume,  we  arrive  at  the  oversimplified  thickness  of  shell  of 
0.22  n  m.  This  is  a  factor  of  two  greater  than  the  more  naive  estimate,  but  it  too  is  only  a  rough 

(albeit  reasonable)  approximation. 

For  the  second  attempt  with  microprobe  equipment,  our  associate,  M.D.  McConnell  of  the 
Microchemical  and  Surface  Analysis  Unit,  Materials  Characterization  Operation,  used  an  alternative 
technique.  The  Ni-Cr  alloy  particles  were  copper-coateu  by  electroless  plating  for  edge  retention, 
and  pressed  together  to  develop  an  apparently  continuous  copper  matrix.  This  compact  could  be 
polished  to  expose  the  particle  cross  section,  and  then  a  traversing  scan  by  Auger  electron 
spectroscopy  might  reveal  an  edge  coating.  The  results  of  such  a  scan  at  8000X  embracing  parts 
of  two  spherical  particles  are  shown  in  Figure  7.  Signals  for  Ni  and  O  are  shown,  that  for  Cr 
being  too  weak  or  overlapping  with  the  others,  There  is  a  clear  buildup  of  oxygen  at  the  particle 
edges,  Its  thickness  could  approach  0.5  /i  m,  but  care  should  be  taken  with  that  interpretation 
without  details  on  the  beam  spot  size. 

The  present  evidence  permits  us  to  conclude  with  confidence  that  there  is  an  oxide  coating 
on  the  particles,  one  which  is  probably  responsible  for  the  high  dc  electrical  resistance  measured 
crudely  on  compacts  of  the  powder  (cf.  next  section).  However,  we  cannot  conclude  the  exact 
chemistry  or  microstructure  of  the  shell  itself,  Conventional  isothermal  oxidation  studies  of  Ni-Cr 
alloys  [15]  show  that  the  nature  of  the  oxide  coating  is  rather  sensitive  to  the  Cr  content  and  to  the 
oxygen  partial  pressure,  as  well  as  to  temperature  and  time.  At  our  level  of  Cr,  it  seems  quite 
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ARBITRARY  UNITS 


possible  that  the  coating  would  contain  precipitates  of  Q2O3  in  a  NiO  host,  but  that  is  still 
conjecture  [16],  The  conditions  of  temperature  and  partial  pressure  of  oxygen  in  the  gas- water 
atomization  spray  are  poorly  known, 


11,5.1  Resistance  Tests  of  Metal  Powder  Compacts 

As  a  final  practical  acceptance  test  for  the  alloy  powder,  we  attempted  to  verify  that  the 
powder  particles  were  electrically  insulated  from  each  other.  Conceptually,  a  check  of  resistance  in 
powder  compacts  appears  simple.  In  practice,  it  is  less  so.  A  simple  device  was  constructed 
consisting  of  solid  copper  electrodes  which  can  slide  into  a  plastic  cylinder,  1.25  cm  inner 
diameter.  We  used  this  with  a  digital  multimeter  capable  of  measuring  up  to  20  MQ  resistance.  A 
small  amount  of  alloy  powder  in  the  cylinder  gave  erratic  results,  dependent  on  pressure,  dme, 
powder  column  height,  etc.  Confirmation  of  this  behavior  was  found  in  the  literature.  Ottavi  et. 
al.  [17]  describe  a  conceptually  elegant  experiment  to  study  percolation  concepts;  i.e.,  at  what 
composition  do  continuous  electrical  paths  develop  when  various  mixtures  of  conducting  and 
insulating  spheres  are  studied?  Their  experiments  were  carried  out  with  monodisperse  5  mm 
diameter  plastic  spheres,  a  fraction  of  which  had  been  electroplated  with  a  thin  copper  coating,  It  is 
helpful  to  quote  from  their  paper,  as  follows:  “In  order  to  obtain  stable  conductance  measurements 
we  find  it  necessary  to  apply  pressure  to  the  piston  bearing  the  upper  electrode.  This  lack  of 
stability  is  not  that  due  to  fluctuations  introduced  by  the  randomness  of  the  mixture.  Indeed,  even 
if  the  system  is  reduced  to  only  two  spheres,  the  measured  resistance  is  not  constant  in  time,  or 
reproducible  between  trials,  unless  a  sufficient  force  (about  10  kgf)  is  applied. . .  We  observe, 
moreover,  a  hysteresis  phenomenon. . ." 

Following  these  hints,  we  obtained  roughly  the  same  behavior.  Powder  column  height 
was  kept  to  1.0  to  1.2  mm.  Significant  arm/body  pressure,  with  some  electrode  piston  rotation, 
was  employed.  These  included  a  hydrogen-reduced  ferromagnetic  metal  powder  and  a  deliberately 
well-coated  magnetic  metallic  powder.  Under  these  conditions,  we  obtained  the  results  given  in 
Table  2.  On  the  basis  of  these  results,  we  conclude  that  the  gas-atomized  Ni-Cr  powder  prepared 
from  the  ingot  Ni92Cr8  should  meet  the  insulated  particle  criterion  for  our  artificial  dielectric. 
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II. 6  Particle  Size  Separations  and  Morphology 

As  mentioned  in  the  Introduction,  particle  size  effects  turn  out  to  be  more  important  than 
originally  anticipated.  Separation  of  -325  mesh  (<37  fj.  m)  Ni-Cr  powder  from  gas  atomization 
preparation  run  RS-67  was  carried  out  (somewhat  tediously)  using  laboratory  sieves  until  about 
200  grams  were  obtained  The  yields  were  as  follows:  29  g  at  <  10  n  m,  65  g  at  10  <  d  <  20  n  m 
and  1 14  g  at  >  20  m.  The  possibility  of  invoking  arch  type  methods  for  millimeter  wavelength 
evaluations  of  composites  at  NRL  necessitated  large  samples  and  hence  greater  quantities  of 
powder.  From  the  later  atomizadon  run  RS-98,  about  1400  g  of  -325  mesh  powder  was  air- 
classified  under  a  subcontract  to  Alpine  American  Corp.  of  Natick,  MA.  The  yields  obtained  were 
as  follows:  550  g  at  <  10  /i  m,  185  g  between  10  and  20  n  m  and  670  g  between  20  and  37  ^  m. 

We  show  some  SEM  photographs  in  Figure  8  of  the  Ni-Cr  powder  from  preparation  run 
RS-67,  for  two  powder  particle  size  ranges,  as  classified  In  each  case  there  are  a  number  of  rather 
smaller  particles  than  called  for  in  the  nominal  range,  but  they  do  not  represent  a  significant  weight 
or  volume  fraction  of  the  batch.  The  other  feature  of  note  is  the  occurrence  of  some  elongated 
particles  among  the  majority  of  spherical  ones.  This  is  more  pronounced  in  the  fraction  with  larger 
sizes.  At  present,  we  have  to  live  with  this,  but  this  fact  is  important  to  know,  both  as  a  reference 
and  as  possible  perturbation  on  the  interpretation  of  data. 


II.7  Preparation  of  Artificial  Dielectric  Composites  -  Polymeric  Binder 
II.7.I  Selection  of  Binder 

One  goal  of  the  program  (Task  1  of  Table  1)  was  the  investigation  of  the  effect  of  binder 
permittivity  in  artificial  dielectrics.  Our  plans  called  for  two  different  insulating  binders,  one  of 
which  would  be  a  polymeric  material  for  which  e  is  typically  3  ±0.5.  The  second  would  be  an 
inorganic  binder  for  which  e  may  be  in  the  range  5  to  10.  For  the  former,  there  is  an  easy  choice 
for  a  laboratory  system,  i.e.,  polyurethane,  known  commercially  as  Permuthane™  (Division  of 
Beatrice  Foods,  Corwin  3t.,  Peabody,  MA.  01960).  It  has  been  used  extensively  in  laboratory 
research  test  composites  by  ourselves  and  our  collaborators  at  GE  Aircraft  Engine,  Evendale, 
Ohio,  who  actually  prepared  our  test  specimen  coupons.  Efforts  at  preparing  a  suitable  inorganic 
binder  are  described  below  in  Section  in. 
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Tabic  2 

RESISTANCE  OF  ALLOY  POWDER  COMPACTS 


Sample-powder 

Approx,  Resistance 

Ni-Cr  <  37/uifl 

I0J-106  n 

Reduced  Magnetic  Metal  <  J7^m 

to1  n 

Well-Coated  Magnetic  Metal  14  <  d  <  27Mm 

>  io7  n 

Table  3 

DESCRIPTION  OF  POLYURETHANE 
TEST  COUPONS  <ALLOY-NlMCr,) 


Sample  # 

Alloy  Particle 
Size 
(in  Mm) 

Nominal  Alloy 
Volume 
Loading,  p 

Measured 
Volume 
Loading,  p 

ONR  1 

- 

0 

— 

ONR  2 

+  20  -  37 

0,40 

0.413 

ONR  3 

+  20  -  37 

0,30 

0.307 

ONR  4 

+  20  •  37 

0,20 

0,200 

ONR  5 

+  20  -  37 

0,10 

0.112 

ONR  6 

+  10-20 

0,40 

0,396 

ONR  7 

-10 

0,40 

0.384 

ONR  8 

-10 

0,20 

0,208 
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II.7.2  Preparation  of  Polyurethane  Test  Specimens 

For  our  initial  series  of  measurements,  a  set  of  nine  test  coupons  was  prepared.  These 
involved  several  particle  sizes,  for  reasons  which  will  be  explained  in  a  later  section.  The  basic 
specifications  are  given  in  Table  3,  including  measured  values  of  p  determined  as  described  below. 

The  preparation  followed  guidelines  for  curing  the  polymer  as  prescribed  by  the 
manufacturer.  Suitably  chosen  amounts  of  alloy  powder  and  binder  powder  were  combined  and 
blended  for  several  hours.  The  mixture  was  leveled  into  a  2.S  cm  x  5  cm  (1  in,  x  2  in,)  die  to  a 
thickness  such  that  the  final  coupon  would  be  -<1  mm  (~0,040  in,)  (In  later  work  somewhat 
thicker  samples  were  used  for  lower  frequency  measurements.)  The  loaded  cold  die  was  placed  on 
hot  press  platens,  then  with  increased  pressure,  raised  to  a  curing  temperature  of  177*C  where  it 
was  held  for  5  min.  It  was  subsequently  cooled  under  cold  running  water. 

The  weighing  calculations  were  based  on  a  cured  Permuthane  density  of  1.06  to  1.1 1 
g/cm3.  (For  sample  ONR  1,  we  determined  1.11  g/cm3.)  These  densities,  as  well  as  those  for  the 
composites  are  obtained  from  dimensions  on  regular  shapes  and  weight.  This  avoids  the  possible 
problems  of  porous  media,  etc.,  but  does  permit  small  uncertainties  in  the  determinations.  For  the 
alloy,  we  used  a  density  of  8.79  g/cm3  based  on  the  magnetic  determination  of  composition  and  an 
x-ray  lattice  parameter  value  of  ao  ■  3.329  ±  0.006A. 


II.7.3  Evaluation  of  Volume  Loading 

As  indicated  in  Section  II.  1.1  we  place  emphasis  on  a  straightforward  and  non-destructive 
approach  to  determining  the  volume  loading,  p ,  of  metal  in  the  artificial  dielectrics,  as  a  guide  to 
analysis  of  their  electromagnetic  properties.  Our  approach  to  this  is  through  magnetic 
measurements  at  low  temperature,  having  chosen  an  alloy  metal  that  has  its  Curie  point  below 
room  temperature. 

A  few  simple  formulas  are  helpful  here.  The  subscripts  a,b,c  and  h  will  refer  to  alloy, 
binder,  composite  and  holes  or  voids  respectively,  Although  we  use  the  symbol  p  to  designate 
volume  loading  in  the  bulk  of  this  report,  for  this  discussion  we  shall  use  the  symbol/with  the 
superscript  w  or  v  to  indicate  weight  fraction  or  volume  fraction.  The  symbols  cr,  M  and  P 
designate  magnetization  per  gram,  magnetization  per  unit  volume  (c.g.,  cm3)  and  specific  gravity 
(or  density  in  g/cm3). 
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I miic  K.  Scanning  electron  nitcrouraphs  of  Nl-Cr  alloy  powder#  from  huh  atoml/.utlon 

preparation  run  RS-67,  IMiiHnlfliutlon;  1000X.  Top  photo,  nomlnul  10  pm  <  d  < 
20  pm;  bottom  photo,  nominal  20  pm  <  ti  <  M  pm. 


2i) 


For  specific  magnetization,  or  magnetization  per  gram,  of  composite  we  write 

but  obviously 

ah  ■  0  and  /"  at  all  times,  and  thus  /* + /,*  ■  1.  Therefore,  we  have 


This  expression  becomes  trivial  but  useful  when  one  notes  that  ak  ■  0  in  most  cases  of  interest. 
Therefore  the  formula  for  weight  fraction  of  alloy  becomes  f?  moJ  cry,  The  measurement  of  a, 
along  with  a  knowledge  (or  independent  measurement)  of  the  alloy  magnetization  gives  its  weight 
fraction. 

For  the  volume  magnetization  and  the  related  volume  fractions,  we  obtain  in  the  same 

manner 


but  this  time  /*”  may  be  different  from  zero  although  Mh »  0  at  all  times  and  Mk  ■  0  in  most  cases 
of  interest  Note  also  that  fi+fi  +/»¥*1,  by  definition.  The  formula  for  volume  fraction  of 
alloy,  p,  is  thus 


p-/: 


ML  m  BtSi.  m  fw  Bl 
Mt  p.o.  h  p. 


This  requires  the  additional  measurement  (or  knowledge)  of  the  density  of  the  composite  and  of  the 
alloy.  With  these  measured  values,  on  whole  samples,  the  results  are  on  firm  ground  and  contain 
no  approximations  or  idealizations.  We  discuss  determination  of  the  void  fraction  In  the  next 
section. 
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As  for  the  case  of  the  binder  alone,  we  have  obtained  more  consistent  results  in  determining 
density  in  the  “primitive"  way  from  dimensions  on  regular  shapes  and  weight.  For  many 
electromagnetic  considerations  it  is  the  overall  volume  property  which  is  of  interest,  so  that  this 
simple  method  which  includes  voids,  cracks,  etc.,  is  the  proper  one.  The  results  for  the  volume 
loading  determinations  are  given  in  Table  3.  The  agreement  between  nominal  and  measured  values 
is  very  satisfactory  and  should  be  characteristic  of  polymeric-type  binders  when  appropriate  care  is 
taken.  With  other  categories  of  binders  such  agreement  cannot  be  taken  for  granted. 


II.7.4  Calculation  of  Porosity 

An  aspect  of  the  microitructure  of  composites  is  the  porosity,  i.e.,  the  volume  fraction  of 
holes  or  voids.  It  is  straightforward  to  calculate  this  quantity  using  the  densities  of  the  separate 
components  and  their  proportions  ascertained  from  observable  measurements.  We  refer  back  to 
definitions  of  symbols  used  above  for  evaluation  of  volume  loading. 

Consider  a  mass  mc  of  the  composite  occupying  volume  vc.  Then 

v,p.  -  mt  •  v,(p./; +pk/; +pj;), 
but  pk  -Oof course.  Thus  we  have 


P, -P.//+PX  or 


■  Pt  ~Pifi 
P* 


But  with  ft  - 1  -  ft  -  ft,  we  can  write 


Pk  Pk  Pk  \Pk  J 


Alternatively,  using  the  formula  for  p  in  observables,  we  obtain 


/;  b1_£l  +  ^A 
Pk  <*. 


,Pk  ”p.. 
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The  principal  sources  of  error  in  this  determination  of  porosity  are  the  density  values  for  the 
binder  material  and  for  the  composite.  By  contrast,  the  alloy  density  and  the  various 
magnetizations  ton  a  mass  basis'!  are  sufficiently  well  known  or  well  measured  to  be  much  less 
important  in  estimating  the  possible  errors  in  .  Therefore,  we  can  write 


+  Qdp, 

dpk  dpt 


fl-.il 

'bJi-sA+ba 

Pk 

*.J 

Pk 

£j£i 

P. 


When  we  substitute  typical  numbers  into  this  formula  we  And:  at  the  p  ■  0.4  level,  a  2%  error 
(uncertainty)  in  contributes  about  0.01  to  and  a  1%  error  in  p,  contributes  a  similar 
uncertainty  (0.01)  to  <#?;  at  the  p  -  0,1  level  each  of  such  a  pair  of  uncertainties  (0.01)  in  dfH 
follows  from  1%  errors  in  or  pt . 

In  Table  4  we  present  the  result  of  the  porosity  calculation.  Other  characteristics  of  these 
samples  such  as  volume  loading  and  particle  size  were  given  in  Table  3.  The  values  of  porosity  are 
generally  small.  The  scatter  between  positive  and  negative  values  (the  latter  are  meaningless)  say 
more  about  our  uncertainties  than  about  the  composite  samples.  At  this  stage  we  conclude  that  the 
preparation  method  is  sufficiently  in  control  that  porosity  is  not  a  significant  factor  for  this  set  of 
samples  to  affect  the  electromagnetic  parameters. 


II.7.5  Follow-on  Higher-Frequency  Test  Specimens 

We  contracted  to  obtain  measurements  at  35  GHz  at  GE  Re-Entry  Systems  Department  in 
Philadelphia.  The  method  was  based  on  the  slotted-line  technique  with  a  circular  waveguide 
cavity.  It  has  the  advantage  of  requiring  a  small  (0.250  in.  diameter),  thin  disk.  This  is  much  less 
demanding  of  material  than,  for  example,  that  required  for  the  arch  method.  It  is  also  easier  to 
fabricate  than  the  precision  (7  mm)  coaxial  toroid,  when  the  question  of  machinability  arises. 

According  to  the  provisions  of  the  contract,  millimeter  wavelength  region  measurements  for 
this  study  are  to  be  performed  at  the  Naval  Research  Laboratory,  Washington,  DC.  Toward  this 
goal,  a  visit  was  made  to  NRL  in  Februaiy  1985  for  discussions  with  Dr.  D.W.  Forester  and  Dr. 
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Tabt*  4 

CALCULATED  POROSITIES 
OF  POLYURETHANE 
TEST  COUPONS 


Sample  # 

rhv 

2 

+  0.04 

3 

-0.03 

4 

-0.03 

5 

+0.01 

6 

-0.04 

7 

+0.07 

8 

+0.01 

Table  S 

DESCRIPTION  OF  POLYURETHANE  TEST  COUPONS 
FOR  NRL  HIGH-FREQUENCY  MEASUREMENT  SAMPLES 


•  Powder  from  atomization  run  RS-98  (Alloy  8.2  «l%  Cr  In  Nl;  7V*  154  K  ±20) 
a  All  Intended  for  nominal  40%  volume  loading  ( p  -  0.40) 
e  Two  powder  particle  size  ranges  chosen 


Coupon  Slab 

Particle  Size 
(Mm) 

Measured 

P 

Measured 

Density 

Calculated 

Porosity 

n 

A 

Waveguide  pieces;  Ka,  Q,  W 

<  10  Mm 

0.36 

3.634 

0.19 

B 

Arch  slab 

<  lOMtn 

0,27 

3.I00 

0.07 

C 

Waveguide  pieces:  K„,  Q,W 

>20-3?Mm 

0.42 

4.31, 

0.01 

D 

Arch  slab 

>  20-37Mm 

0,39 

4.040 

0.08 
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FJ.  Rachford.  It  was  agreed  that  the  measurements  would  be  attempted  by  several  techniques, 
including  the  arch  method  at  rather  high  frequencies  (~90  GHz),  and  waveguide  techniques  at 
several  high-frequency  bands,  e.g.,  Ku  (26-40  GHz),  Q  (40-60  GHz),  and  W  (75-1 10  GHz). 
Such  a  wide  range  of  frequencies,  coupled  with  the  particle  size  effects,  offers  an  exciting  prospect 
for  testing  our  models  and  amplifying  our  understanding. 

From  the  discussions  at  NRL,  it  was  decided  to  measure  rather  highly  loaded  samples,  i.e., 
volume  loading  of  p  -  0.40.  In  order  to  preserve  the  pardcle  size  study,  we  chose  the  two 
extremes,  i.e.,  <  10  n  m  and  20  to  37  ji  m.  Test  coupons  of  powder  (RS-98)  with  polyurethane 
binder  were  prepared  according  to  the  above  techniques,  Two  coupons  for  each  particle  size  were 
made.  One  was  4  x  4  x  0.045  in.  for  preparation  of  the  arch  method  test  piece,  and  the  second  was 
about  1  x  2  x  0.080  in.  for  subsequent  precision  cutting  of  waveguide  pieces. 

The  actual  volume  loading  and  estimated  porosities  were  calculated  as  described  above. 

The  results  are  presented  in  Table  5.  We  see  that  the  larger  powder  size  material  coupons  came 
quite  close  to  the  goals  set  (p  *  0.4)  while  the  <  10  n  m  powder  composite  coupons  fell  short. 

At  millimeter  wavelengths,  the  samples  must  be  cut  with  high  precision.  The  3  x  3-in. 
plate  samples  for  the  arch  method  were  trimmed  to  thicknesses  estimated  at  (3/4)  Xm ,  where  Xm  is 
the  wavelength  in  the  medium  at  W  band  (90  GHz).  Estimates  for  Xm  were  arrived  at  from 
modeling  for  /i'  and  e'  described  herein.  The  formula  for  the  wavelength  in  the  medium  is  Xm  - 

Xt  /  [eV/]1/l*  where  Xt  is  the  vacuum  wavelength.  Even  more  demanding  are  the  samples  to  fit 
into  the  waveguides  of  the  separate  bands  These  desirable  thicknesses  are  again  an  odd  number  of 
quarter- wavelengths.  However,  for  this  measurement  one  refers  to  the  guide  wavelength  [4],  Xt 
given  by 


X 


i 


where  Xe  is  the  cutoff  wavelength  of  the  particular  waveguide  which  is  a  function  of  the 
waveguide  dimensions  [20].  In  some  cases  we  could  get  by  with  quarter-wavelength  pieces,  At 
higher  frequencies,  these  would  be  so  thin  as  to  force  a  choice  to  three-quarter-wavelength 
thicknesses. 
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II.8  Direct  Observations  on  Composites:  Porosity  and  Particle  Aggregation 

Can  we  confirm  the  calculated  conclusion  that  porosity  in  the  polyurethane  artificial 
dielectric  composites  is  negligible?  That  would  seem  to  be  an  easy  task  in  these  days  when  a 
scanning  electron  microscope  (SEM)  is  often  readily  available.  However,  preparation  of  samples 
for  observation  is  a  special  an  which  is  needed  in  the  case  of  these  polymeric  composites.  Special 
equipment  was  needed,  as  provided  by  our  associates  V.H.  Watkins  and  S.Y.  Hobbs  of  the 
Polymer  Physics  Unit.  Room  temperature  microtoming  was  unsatisfactory  because  the  binder  was 
so  soft,  but  a  cold  microtoming  technique  (with  the  sample  at  about  -150°C)  using  a  glass  knife 
produced  a  very  satisfactory  surface  for  examination.  The  apparatus  used  was  an  Ultracut  E. 

For  the  first  sample  studied,  the  as-cut  material  was  lightly  coated  with  Au  and  Pd  (a  few 
hundred  angstroms  coat)  to  carry  off  electronic  charge  and  heat  The  SEM  pictures  shown  in 
Figure  9  for  one  sample  (ONR  #2,  p  «  0,4, 20  /i  m  <  d  <  37  /im)  are  clean  and  interesting,  There 
are  faint  but  discemable  circles  which  are  cut(l)  sections  of  the  metal  powder  spheres.  Their 
dimensions  correspond  to  particle  sizes  used.  Note  that  smaller  sizes  may  be  anticipated  in  the  cut, 
and  a  few  are  distorted  from  the  circular  cross  section.  Their  identity  is  confirmed  by  EDAX 
(Energy  Dispersive  Analysis  by  X-rays)  which  shows  a  high  count  for  Ni  and  a  lesser  one  for  Cr 
on  these  circles,  in  contrast  to  the  adjoining  matrix.  A  number  of  holes  are  present.  Some  of  these 
are  obviously  hemispherical  “pull-outs."  Others  are  distorted  but  with  excess  lip  material 
suggesting  a  pull-out.  There  are  many  regions  densely  filled  with  material  of  one  kind  or  another. 
The  latter  observation  is  in  accord  with  the  calculations  suggesting  negligible  porosity.  At  this 
volume  loading  the  metal  particle  density  is  sufficiently  high  that  one  cannot  get  a  feel  for  “irregular 
fluctuations"  in  the  loading. 

It  is  a  matter  of  concern  for  understanding  the  electromagnetic  properties  on  a  quantitative 
basis  that  one  look  for  wide  variations  in  local  packing  on  some  moderate  length  scale.  Several 
workers  have  pointed  out  the  potentially  enormous  effects  that  may  arise  from  agglomeration  with 
electric  conduction  contact  {18]  or  wide  variation  in  size  distribution  [19]  intrinsically  or  as  a 
consequence  of  agglomeration.  In  addition,  as  we  discuss  later  herein,  if  a  nominally  dilute,  low 
p- value  composite  has  locally  densely  packed  regions  (of  a  size  we  hope  to  be  able  to  investigate 
theoretically),  compensated  by  regions  nearly  empty  of  metallic  particles,  without  electric 
conduction  paths  in  the  densely  packed  regions,  there  can  be  some  enhancement  of  electromagnetic 
property  parameters.  This  arises  from  a  strongly  non-linear  volume  loading  dependence  of  the 
property. 
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re  9.  Scanning  electron  microscope  views  of  cold-microtomed  polymeric  composite 
(lightly  coated  with  Au  and  Pd):  ONR  #2  (p  =  0.4,  20  ttm  <  <1  <  37  nm).  Top 
photo  at  1S0X;  lower  photo  at  500X. 
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Subsequently  we  extended  the  examination  to  the  mote  dilute  samples  of  the  series  (ONR 
No.  5,  p  =  0. 1 ;  ONR  No,  4,  p  *  0.2;  ONR  No.  3,  p  =  0.3;  20  fi  m  <  d  <  37  fJ.  m,  alloy  RS-67). 
This  time  the  microtomed  samples  were  coated  lightly  with  carbon  to  carry  off  electrons  and  heat, 
instead  of  Au  and  Pd.  The  result  is  a  much  improved  visual  contrast,  as  shown  in  Figure  10. 
There  is  clear  evidence  for  spatial  fluctuations  in  local  density  of  alloy  particles.  This  picture  gives 
a  scale  to  the  “clumping.”  The  challenge  is  whether  we  can  model  this  in  an  analytic  way.  The 
SEM  pictures  in  Figure  10a  also  show  “crevasse”  lines  in  the  matrix  that  do  not  appear  in  the  other 
photos,  We  (including  Watkins  and  Hobbs)  do  not  know  their  origin,  but  think  they  are  not 
relevant  to  our  problem.  Our  collaboration  with  Professor  Doyle  was  initiated  in  hopes  of  tackling 
the  clumping  in  a  quantitative  way. 
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Figure  10b.  Scanning  electron  microscope  views  of  cold-mlcrotomed  polymeric  composite 
ONR  #4  (p  -  0.21i  same  d),  Top  photo  at  1  SOX j  lower  photo  at  500X. 
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Figure  10c.  Scanning  electron  microscope  views  of  cold-mlcrotomed  polymeric  composite, 
ONR  #3  (p  =  0.3.  same  </).  Top  photo  at  1  SOX ;  lower  photo  at  500X. 
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Ill,  BINDER  PERMITTIVITY  EFFECTS 

An  aspect  of  the  modeling  analyses  of  the  permittivity  of  artificial  dielectrics  is  that  the 
dielectric  constant  of  the  binder  “sets  the  level”  for  the  dielectric  constant  of  the  composite,  and  in  a 
multiplicative  way,  This  can  be  expressed  quite  simply  by  citing  the  formula  for  the  Mitoff 
approximation  to  the  Maxwell  (Clausius-Mossotti)  formula  [21],  i.e.: 

logs'  ai  log  fij  +  Gp 


or 


£,'  /  e'b  a  10  Op, 

where  c  and  b  represent  composite  and  binder,  p  is  the  volume  loading  of  the  metal  or  alloy 
particles  (phase  a),  and  G  is  a  constant  of  value  1.3  Titus  choosing  a  binder  with  twice  or  half  the 
value  for  will  produce  a  doubling  or  halving  the  value  of  e'. 

The  relevance  of  this  effect  lies  in  the  role  of  £«  in  the  wavelength  of  the  electromagnetic 
radiation  in  the  medium,  by  the  formula  given  in  Section  D.7,5.  Since  electromagnetic  applications 
configurations  often  are  designed  with  materials  of  thickness  measured  in  units  of  the  wavelength, 
the  choice  of  binder  could  influence  the  overall  weight  of  the  configuration. 

As  noted  in  Section  n.7.1,  one  of  the  tasks  was  to  verify  this  prediction  in  the  microwave 
and/or  millimeter  wavelength  regions  by  comparing  artificial  dielectric  composites  with  two 
different  insulating  binders.  (In  hindsight,  the  proposed  task  was  rather  naive  in  that  the  prediction 
has  been  rather  well  established.)  Our  plan  called  for  one  organic,  polymeric  binder,  for  which  £ J 
is  usually  near  3,  and  a  second  inorganic  binder  that  could  carry  an  e't  value  in  the  range  of  5  to 
10. 


IIM  Criteria 

At  the  outset  it  is  useful  to  enumerate  the  requirements  of  the  proposed  search.  The 
material  and  process  should  maintain  the  shape  (usually  spherical)  of  the  alloy  particles;  they 
should  maintain  the  essential  chemical  integrity  of  the  alloy,  although  some  small  deviation  is 
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permitted;  and  they  should  maintain  a  condition  of  electrical  isoladon  of  the  particles.  Beyond 
these  requirements,  the  resulting  product  should  have  sufficient  mechanical  integrity  to  permit 
machining  into  high-frequency  measurement  specimens.  Even  with  the  last  requirement,  there  are 
options  in  that  precision  coaxial  samples  for  the  microwave  network  analyzer  are  more  demanding 
of  machinability  than  the  thin,  solid  disks  needed  for  the  35-GHz,  slotted-line  measurement. 


in.2  Extensive  Search  for  Alternative  Inorganic  Binder 

For  our  first  attempt  we  tried  to  use  a  pioprietary  inorganic  binder  under  development  at 
GB-Evendale.  It  is  intended  for  a  similarly  sized  but  chemically  different  dispersoid  (loading 
material).  Typically,  with  inorganic  binders  there  is  a  high-temperature  cycle  to  bring  the  binder  to 
its  final  state.  In  our  past  experience  [1]  such  a  high-temperature  excursion  may  be  accompanied 
by  a  chemical  reaction  between  binder  and  dispersoid  particle.  In  the  present  case  we  requested 
four  composites,  p  »  0.1, 0.2, 0.3,  and  0.4,  each  with  the  largest  particle  size  +20-37  p  m  of  the 
Ni92  Cub  powder,  In  each  case,  the  resulting  composite  was  ferromagnetic  at  room  temperature. 
Indeed  a  thermomagnedc  scan  showed  Curie  points  up  to  350*C  (pure  Ni,  Tc  *  358'C)  with  strong 
suggestions  of  composition  variations  within  the  magnetic  phase.  Under  these  conditions  it  is  very 
difficult  to  know  the  alloy  composition  in  the  composite.  Therefore,  we  could  not  determine 
volume  loading  without  destruction  and  visual  or  probe  analysis  of  the  composite. 

The  subsequent  steps  of  our  search  were  carried  out  by  our  collaborators  M.B.  Borom  and 
L.E.  Szala.  Specimens  were  compounded  from  a  three-component  system-paramagnetic  Ni-Cr 
alloy  particles  as  the  active  component,  a  choice  of  two  different  glasses  as  the  binder,  and  alumina 
as  an  inert  filler.  Ni-Cr  alloy  particles  of  two  different  particle  sizes  were  used  in  volume  fractions 
ranging  from  13  to  30  vol%.  Alumina  as  an  inert,  refractory,  fine-grained,  dielectric  material  was 
chosen  to  keep  the  metallic  particles  isolated  from  one  another,  Good,  high-temperature 
compressive  strength  of  the  alumina  contributes  to  maintaining  spacing  of  the  deformable  metallic 
phase,  particularly  during  sintering  or  hot  pressing.  The  third  component  to  be  considered  is  the 
binder.  (Electrically,  all  the  nonmetallic  components  make  up  the  binder,  but  we  are  speaking  here 
in  the  ceramic  processing  sense).  The  binder  must  be  sufficiently  fluid  during  processing  to  flow 
between  the  AI2O3  and  the  alloy  and  create  a  mechanically  strong  structure.  The  binder  must  flow 
at  a  temperature  below  the  melting  point  of  the  Ni-Cr  alloy,  and  must  also  be  nonreactive  to  it. 
Many  glasses  with  reasonable  softening  points  fit  these  requirements.  Two  types  of  glass  were 
tried:  a  soda/lime  window  glass  composition  designated  as  Coming  code  0800,  and  a  lead- 
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bearing,  Owens-Illinois  solder  glass  designated  as  SG-7.  Both  glasses  have  softening  points 
below  the  melting  point  of  the  Ni-Cr  alloy. 

Powder  preparation  was  accomplished  by  milling  the  components  in  polyethylene  jars 
using  1/4-inch  alumina  spheres  as  the  grinding  material.  The  alumina,  0.3  Jim  in  size,  and  the 
glass  were  initially  wet-milled  together  in  alcohol  for  8  to  10  hours  to  reduce  the  panicle  size  of  the 
glass  and  to  ensure  proper  blending  of  the  two  components.  After  the  addition  of  the  alloy 
powder,  an  additional  1  to  2  hours  of  milling  was  performed  to  disperse  the  alloy  powder.  Milling 
times  after  the  addition  of  alloy  were  kept  shon  to  prevent  erosion  of  the  oxide  coating  already 
present  on  the  alloy  powder  particles.  The  blended  powders  were  then  dried  in  a  rotating 
evaporator  to  keep  the  relatively  coarse-grained  Ni-Cr  in  suspension.  The  dried  blends  were 
compacted  using  standard  die  and/or  isostatic  pressing. 

Table  6  is  a  summary  of  experiments  and  results.  The  table  is  made  up  of  two  sections, 
with  the  type  of  glass  used  being  the  prominent  variable.  Samples  are  listed  alphabetically  in  each 
section.  The  m^jor  sections  are  then  divided  into  segments  detailing  composition  and  process 
variables  and  results. 

Sample  0800-A,  a  blend  of  35%  AI2O3, 35%  0800  soda/lime  glass  and  30%  20-37  n  m 
Ni-Cr  powder,  was  consolidated  by  hot  pressing  at  a  temperature  of  950*C  and  600  psi  loading, 
The  resultant  material  shown  in  Figure  1 1  appears  to  meet  the  requirements  of  paramagnetism  (i.e„ 
not  attracted  by  a  small  permanent  magnet)  and  is  nonconductive;  however,  the  compact  has  very 
low  strength. 

In  an  attempt  to  overcome  the  low  strength  problems  associated  with  sample  0800-A,  a 
higher  hot  pressing  temperature  and  loading  were  employed.  It  was  felt  that  increasing  the 
temperature  to  1300'C  and  the  pressure  to  1200  psi  would  require  a  more  rigid,  less  deformable 
composition.  This  Ni-Cr  content  was  kept  at  30  vol%,  while  the  glass  content  was  decreased  to 
only  10%,  and  the  difference  made  up  by  alumina  (sample  0800-B),  As  anticipated,  use  of  higher 
temperatures  and  pressures  corrected  the  low  strength  problem  associated  with  sample  0800-A,  but 
created  a  problem  of  its  own.  Because  of  the  elevated  temperature  and  pressure,  the  alloy  became 
more  deformable  and  sintered  into  an  interconnected  network,  thus  becoming  electrically 
conductive. 

Sample  0800-C  was  prepared  keeping  the  ratio  of  AI2O3  to  glass  the  same  as  in  0800-B, 
but  reducing  the  metallic  alloy  loading  to  one  half  the  value,  i.e„  15%.  Hot-pressing  conditions 
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Klpure  11.  Optical  mlcrojiraphs  of  lnor|(nnic  composite  nii  in  pie  OHOO-A,  with  soda-lime 
pluss,  Preparation  parameters  given  In  Table  6.  Magnifications  at  200\  and 


were  identical  to  0800-B.  These  conditions  produced  a  very  dense,  paramagnetic,  nonconductive 
sample.  The  process,  however,  would  probably  not  be  suitable  because  of  the  low  loading  of  Ni- 
Cr  alloy. 

Sample  0800-D  differs  from  sample  0800-B  only  in  that  the  alloy  particle  size  was  reduced 
from  20-37  to  <10  fj.  m.  Composition  and  processing  were  the  same.  A  smaller  particle  size  was 
used  in  an  attempt  to  increase  the  loading  back  to  30%  and  to  achieve  a  better  dispersion  of  the  Ni- 
Cr  alloy  throughout  the  mixture,  thus  possibly  creating  less  contact  between  the  individual 
particles.  It  is  very  difficult  to  suspend  large  particles  of  a  dense  material  in  a  fine  dispersion  of  a 
low-density  material.  Results  of  the  test  were  poor,  Even  though  the  uniformity  of  the  dispersion 
of  Ni-Cr  alloy  was  greatly  increased,  the  material  was  highly  conductive.  The  result  is  shown  in 
Figure  12 . 

The  Ni-Cr  alloy  powder  appears  to  be  too  deformable  at  elevated  temperatures  to  withstand 
unidirectional  pressure  of  hot  pressing.  Pressure  exerted  on  the  system  deforms  the  spherical 
metal  particles  into  platelets  and  increases  the  chance  of  contacting  and  sintering  together. 
Therefore,  experiments  employing  pressureless  sintering  were  designed  in  an  attempt  to  solve  the 
problem  of  metal-metal  contact  during  consolidation. 

Samples  0800-E  and  0800-F  were  compositionally  the  same  as  0800-D  (60%  AI2O3, 30% 
Ni-Cr  (<10  n  m),  10%  0800  glass).  The  mixtures  were  isostatically  pressed  into  1/2-in.  diameter 
x  3/8-in.  long  pellets  at  a  pressure  of  S3  kpsi.  In  an  effort  to  increase  the  oxide  film  thickness  on 
the  Ni-Cr  particles  to  assist  in  preventing  metal-metal  contact,  sample  0800-E  was  subjected  to  a 
preoxidation  step  by  firing  the  sample  to  800'C  in  air  and  holding  for  10  min.  Both  samples  0800- 
E  and  0800-F  were  sintered  by  heating  in  helium  at  ~50'C/min  to  1000'C  and  held  at  temperature 
for  3  min.  Even  though  both  samples  remained  electrically  nonconductive,  the  paramagnetic 
properties  of  each  sample  were  destroyed,  The  loss  of  paramagnetism  and  development  of 
ferromagnetism  were  indicated  by  the  ease  with  which  the  sample  could  be  picked  up  and  held  by  a 
small  permanent  magnet,  Loss  of  paramagnetic  properties  was  not  observed  in  any  of  the  hot- 
pressed  samples,  The  change  in  magnetic  properties  indicates  that  the  Ni-Cr  alloy  had  been 
changed  compositionally,  most  likely  due  to  oxidation  during  sintering  and/or  preoxidation. 

From  both  sintering  and  hot-pressing  experiments  it  became  apparent  that  temperatures  in 
excess  of-900'C  could  be  detrimental  to  the  Ni-Cr  alloy  particles  both  chemically  as  well  as 
physically,  To  reduce  the  processing  temperature  below  900‘C,  a  lower  softening  point  glass  was 
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glasn.  Preparation  parameters  given  in  Table  6.  Magnifications  at  200X  and 
500X, 
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selected.  The  glass  selected  was  a  lead-bearing  solder  glass  designated  SG-7,  which  has  a 
working  temperature  of  approximately  600-700’C. 

Sample  SG-7  A,  a  mixture  of  79%  AI2O3  and  21%  SG-7  containing  no  metal  particles,  was 
hot  pressed  under  a  static  load  of  700  psi.  The  sample  completely  densified  at  690*C. 

Sample  SG-7B  was  formulated  with  30%  Ni-Cr  (20*37  fi  m).  The  volume  percent  of  SG- 
7  was  kept  at  the  same  level  as  in  SG-7A  (20%)  and  the  balance  was  alumina.  The  mixture  was 
hot-pressed  under  a  static  load  of  400  psi.  In  contrast  to  sample  SG-7  A,  which  glassifled  at 
690*C,  no  movement  of  the  press  rams  occurred  until  750'C.  The  final  pressing  temperature  was 
840'C.  Both  the  reduced  pressure  on  the  sample  and  the  presence  of  the  coarse  metal  particles 
probably  contributed  to  the  need  for  higher  processing  temperature.  The  density  of  the  resultant 
compact  was  excellent.  The  material  was  non-conductive  and  maintained  its  paramagnetic 
properties,  but  as  in  hot-press  experiments  the  alloy  was  no  longer  spherical. 

With  further  evidence  that  hot-pressing,  even  at  850‘C,  deforms  the  Ni-Cr  alloy  particles, 
sintering  was  again  employed.  Sample  SG-7C  is  compositionally  identical  to  SG-7B.  The 
material  was  sintered  at  a  temperature  of  1000'C  for  3  min  in  1  atm  of  helium.  The  high 
temperature  was  again  tried,  even  though  previously  unsuccessful,  as  a  direct  comparison  between 
the  0800  and  the  SG-7  glass  (via  sample  0800-E).  The  results  were  nearly  the  same.  The 
paramagnetic  properties  were  destroyed.  The  material  is  nonconductive,  but  this  is  primarily 
because  the  material  is  substantially  unsintered.  A  higher  content  of  glass  may  be  required  for 
densification. 

After  discussions  it  was  decided  to  relax  the  potential  high-temperature  service 
requirements  of  the  material  to  see  if  Ni-Cr  could  be  successfully  incorporated  in  a  low-temperature 
glass  matrix  without  the  presence  of  the  refractory  alumina  phase. 

Sample  SG-7D  is  a  mixture  of  70%  SG-7  solder  glass  and  30%  Ni-Cr  alloy  (20-37  m). 
The  intent  of  the  experiment  was  to  sinter  at  -  630*  C,  but  due  to  an  error  in  reading  a  thermocouple 
scale,  the  temperature  was  raised  to ~ 1000'C,  Again,  and  not  surprisingly,  the  sample  failed  all 
the  necessary  requirements  (i.e.,  it  was  both  ferromagnetic  and  electrically  conductive). 
Micrographs  are  shown  in  Figure  13. 

The  experiment  was  then  repeated,  this  time  sintering  at  625’C  for  20  min.  The  results 
(samples  SG-7E)  were  quite  encouraging.  The  paramagnetic  and  electrical  conductivity 
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13.  Micrographs  of  Inorganic  composite  sample  S(»7-l),  with  lead-solder  glass. 
Preparation  parameters  given  in  Table  6.  Magnifications  at  200X  and  500X. 


requirements  appear  to  have  been  met.  Microstructural  observations  in  Figure  14  show  the 
material  to  be  quite  dense.  The  Ni*Cr  powder  has  remained  substantially  spherical,  when 
compared  with  its  starting  state  in  Figure  8. 

III.3  Initial  Evaluation  Measurements 

The  procedures  used  to  make  sample  SG-7E  were  continued  so  as  to  prepare  a  series  of 
three  sample  pieces  with  a  range  of  nominal  volume  loading  (p)  values.  These  were  evaluated  first 
by  thermomagnetic  analysis,  similar  to  that  shown  in  Figure  2,  which  first  verified  that  the  desired 
characteristic  of  no  ferromagnetism  at  room  temperature  was  retained,  and  second,  provided  data 
on  Curie-temperature  and  low-temperature  (6K)  saturation  magnetization  from  which  the  actual 
volume  loading  could  be  calculated,  The  results  are  presented  in  Table  7.  The  Curie  temperatures 
allow  one  to  determine  the  alloy  composition  from  Figure  3.  The  saturation  magnetization  of  the 
Ni-Cr  alloy  present,  in  the  composite  follows  from  Figure  4,  and  when  combined  with  densities  of 
alloy  and  composite  the  volume  loading  is  obtained  from  p  =  aepe  /  crjv 

The  resulting  values  are  quite  acceptable  for  the  nominal  20%  and  30%  loadings.  By 
contrast,  that  for  the  nominal  40%  samples  is  poor.  Its  low  density  correlates  with  a  microscopic 
observation  of  significant  porosity.  Another  interesting  point  appears  in  Table  7:  The  Curie 
temperatures  are  all  slightly  lower,  by  15  to  20K,  than  those  shown  in  Figure  2  for  the  original  as- 
atomized  powder.  If  we  assume  that  there  was  no  sampling  problem,  nor  a  significant  variation  of 
Tc  with  particle  size,  we  are  led  to  conclude  that  a  slight  reaction  did  occur  in  preparing  the 
composite,  which  shifted  the  Ni-Cr  ratio  slightly  toward  that  of  the  original  ingot.  This  would  not 
impair  the  application  of  this  material  for  our  purposes, 

The  other  “evaluation”  technique  attempted  was  the  preparation,  by  careful  machining,  of 
precision  coaxial  toroids  for  high  frequency  measurements.  The  results  were  not  a  success,  in  that 
the  composites  were  very  brittle  and  tended  to  chip  or  shatter.  Where  two  toroids  were  needed  for 
measurement  at  each  volume  loading,  we  were  lucky  to  obtain  one,  and  that  one  might  have  had 
small  chips  at  its  corners.  Therefore,  we  confined  this  series  to  the  35-GHz  measurements,  which 
require  only  a  single  solid  disk  for  the  slotted-line  technique. 


Figure  14.  Micrographs  of  inorganic  composite  sample  SG7-E,  with  ieud-solder  glass. 

Preparation  parameter:,  given  In  Table  6.  Magnifications  nt  100X  and  500 X. 


I 
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Table  7 


THERMOMAGNETIC  EVALUATION  PARAMETERS 
OF  INORGANIC  COMPOSITES  PREPARE*)  WITH  LEAD-SOLDER  GLASS 
BY  LOW  TEMPERATURE  SINTERING 


Semple  Idem 
(Powder  Source) 

Nominal 

Alloy 

(vol%) 

Curie 

Temp, 

Sat.  Mag. 

<fm  (6  K) 

(emu/g) 

At% 

Cr 

Specific 

Gravity 

Measured 
Vol  Load 
(P) 

SG-7-S 

(RS-67;  10<4<20) 

20 

143  K 
±21 

5.30 

8.2 

4.47 

0.18 

SG-7-E 

(RS-98;  20<4<37) 

30 

137  K 
±23 

6.75 

8.4 

5.01 

0.27 

SG-7-H 

(RS-67;  10<4<20) 

40 

147  K 
±26 

8,97 

8.2 

4.23 

0.28 
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IV.  THEORETICAL  APPROACHES/MODELING 


In  this  section  we  address  the  challenges  posed  in  Task  3  of  the  Statement  of  Work  (Table 
1).  Despite  the  "classical”  background  of  artificial  dielectric  composites  [2]  in  the  nineteenth 
century,  the  subject  remains  quite  alive  today.  There  have  been  several  conferences  in  recent  years 
on  the  properties  of  heterogeneous  media  [22,23].  In  addition  to  the  obvious  dielectric  or 
permittivity  behavior,  this  particular  class  of  mixture  also  can  exhibit  permeability  phenomena.  We 
begin  our  discussion  on  this  point, 


IV, 1  Magnetic  Dipole  and  Particle  Size  Effects 

At  first  consideration  it  might  seem  that  purely  electric  dipole  phenomena  would  fully 
describe  the  high-frequency  response  of  an  artificial  dielectric  composed  of  strictly  non-magnetic 
materials  (in  contrast  to  prior  work  on  artificial  magneto-dielectrics  [1]).  Indeed  this  appears  to 
have  been  the  initial  point  of  departure  for  analysis  over  the  past  rather  active  decade  of  research  on 
selective  absorbers  of  solar  and  infrared  radiation  [24, 22].  However,  it  was  quickly  realized  that, 
if  the  diameter  of  the  conducting  particle  is  not  negligible  compared  to  the  skin  depth,  there  is  a 
second,  rather  classical,  mechanism  for  absorption  arising  from  the  effects  of  the  magnetic  field  of 
the  incident  radiation.  In  the  time-varying  magnetic  field,  there  will  be  an  induced  eddy  current 
magnedc  polarization  and  an  eddy  current  loss.  The  effective  permeability  of  the  composite 
artificial  dielectric  will  thus  not  be  unity  even  though  the  dc  permeability  of  the  constituents  may  be 
exactly  unity.  These  magnedc  effects  are  size-dependent  while  the  electric  dipole  polaiizadon  and 
loss  phenomena  are  not.  Indeed  the  magnedc  contribution  dominates  for  metallic  panicles  i*bove 
about  30  A  in  size  [24, 25]! 

As  noted  in  Section  1,  this  feature  was  not  in  our  original  plans.  The  calculadons  presented 
below  show  that  much  interesting  behavior  occurs  in  relevant  frequency  and  panicle  size  ranges. 
Therefore  we  put  emphasis  on  obtaining  and  measuring  different  powder  size  fabrications  (Section 
II).  This  is  the  expansion  of  scope  noted  in  Section  I. 


IV.1.1  Magnetic  Polarizability  of  a  Single  Particle 

The  magnetic  response  of  an  isolated  isotropic  conducting  sphere  has  been  discussed  nicely 
by  Landau  and  Lifshitz  [26].  Following  their  notation, 
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A imtxH 


where  M  is  the  magnetization  (magnetic  moment  per  unit  volume),  H  is  the  external  magnetic  field 
and  a  is  the  (complex)  magnetic  polarizability.  For  a  sphere  they  derive  (with  a  =  a'  +  ia"), 

36[slnh(2fl  /  6)-  sin(2a  /  £)] ' 

2a[cosh(2fl  /  8)-  cos(2o  /  <5)J 


and 


96 8  {  asinh(2a  /  g) + sin(2a  /  6) 
1 6ffa2  _  5cosh(2a/5)-cos(2fl/5) 


where  6  is  the  skin  depth  given  by 


£ 

6  *  iii  i  in  cm  (with  cgs  units) 

^2  ;p(raa> 


or 


<5  . . . . .  in  m  (for  practical  or  SI  units) 

In  these  expressions  cru  ■  1  /  pu  is  the  electrical  conductivity  (in  appropriate  units), /the 
frequency  and  p,  the  permeability  of  free  space  (4tt  x  10*7  henry/m). 

In  Figure  IS  we  plot  -a',  the  (negative)  real  part  of  the  conducting  sphere  magnetic 
polarizability  as  a  function  of  a  /  6,  the  radius  divided  by  the  skin  depth.  In  the  low  frequency  (or 
small  particle  size)  limit,  6  »  a ,  this  component  decreases  as  (a  /  <5)4,  At  the  other  extreme,  it 
reaches  a  finite  value,  af  *  -3/8 ff,  corresponding  to  the  magnetic  behavior  of  a  superconducting 
sphere  [26]. 
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(MAGMARY  RWT)«“ 


(o/l)  RADIUS  /SKIN  DEPTH 

Figure  IS.  Magnetic  polarliablllty  of  a  conducting  aphere.  Negative  real  part,  -a',  vi.  the 
dlmenilonleu  ratio  of  sphere  radius  to  skin  depth. 
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Figure  16.  Magnetic  polarliablllty  (imaginary  or  loss  part)  of  conducting  sphere  In  time 
varying  field,  a"  vs.  {a  16) 
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In  Figure  16,  the  imaginary  part  of  the  magnetic  polarizability  is  presented,  again  as  a 
function  of  a  /  S.  This  quantity  is  intrinsically  positive.  It  vanishes  in  both  of  the  extremes  of 
small  or  large  radius  (or  frequency);  in  the  former  the  limiting  behavior  is  as  (a  /  <5)2,  while  in  the 
latter  it  decreases  in  proportion  to  (S  /  a).  The  maximum  loss  occurs  at  a  /  S  -  2.3,  while 
remaining  above  half  its  maximum  for  the  (a  /  S)  range  of  1.2  to  7. 


IV.  1.2  From  Single  Conducting  Sphere  to  Composite 

The  magnetic  polarizability  formulas  given  in  the  preceding  section  pertain  only  to  the 
spherical  shape  and  the  externally  applied  magnetic  field.  In  describing  magnedc  properties  in 
general,  as  well  as  in  the  present  case  of  a  dispersion  of  small  magnetic  particles,  one  must  proceed 
carefblly  in  accounting  for  the  various  demagnetizing  effects,  internal  and  external  fields,  etc.  For 
example,  in  the  present  case  one  may  have  to  distinguish  between  the  demagnetizing  coefficients  of 
the  overall  composite  sample,  of  the  individual  particle,  and  of  the  cavity  “constructed"  around  the 
particle  in  analyzing  for  interactions  between  the  particle's  dipole  moment  and  the  magnetization  of 
the  particles.  In  our  model,  the  latter  interactions  are  considered  in  the  mean-field  approximation, 
which  is  accurate  in  the  dilute  limit.  (We  shall  describe  some  of  its  weaknesses  in  a  later  section 
when  we  deal  with  the  dielectric  behavior).  As  another  example,  particularly  relevant  to 
conventional  magnetic  materials,  i.e.,  ferromagnetic  ones,  the  quantity  of  fundamental  importance 
is  the  intrinsic  or  internal  property,  be  it  susceptibility,  permeability,  or  polarizability. 

For  the  artificial  dielectric  composites  studied  here,  we  have  the  simple  case  of  spherical 
particles  (demagnetizing  coefficient  ■  4*3 ),  spherical  cavity  (same),  and  overall  sample 
configurations,  for  which  we  can  ignore  demagnetization,  The  latter  come  about  because  the 
measurements  are  carried  out  with  closed  flux  configurations  (toroid  in  a  co-axial  line)  or  thin  slabs 
with  the  incident  H-fleld  parallel  to  the  slab  plane. 

Having  hoisted  the  flags  of  caution,  we  shall  not  repeat  the  derivations  which  may  be  found 
in  appropriate  textbooks,  but  shall  merely  present  some  key  formulas.  The  intrinsic  (or  internal) 
magnetic  susceptibility  of  the  sphere  (relative  to  the  internal  Held)  is  related  to  the  polarizability 
derived  above  as; 


1  -  4/ra  /  3 
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For  a  composite  of  spherical  particles  with  volume  loading,  p,  the  susceptibility  relative  to  an 
externally  applied  field  is: 


tail 


_ XiE _ _ 


where  Ns  is  the  demagnetizing  coefficient  for  the  overall  composite  sample  shape.  In  the  common 
case  herein,  we  have  Nt  ■  0,  and  thus  we  obtain 

y  -  .  PXl 

“  ufe-pj*  ' 

Upon  substituting  the  relation  between  &  and  a,  this  formula  becomes 

't* 


This  appears  so  simple  that  one  presumes  it  could  have  been  found  more  directly.  However,  we 
feel  the  excursion  into  details  and  warnings  is  appropriate.  Finally  for  the  permeability,  we  have, 
by  definition,  *  1  +  and  then 


At  this  point  we  take  note  that  a  and  hence  p  are  complex,  with  the  forms  a  *  o' + lot*, . 
fi  ■  p'  +  iji*.  Making  these  substitutions  and  continuing  the  manipulation,  we  obtain 


1  + 


4?rpg* 

D 
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with 


'  D-l-2(lifjpa'+(^-pj[{af +(«")•]  . 

In  Figures  17  and  18  we  show  the  behavior  of  M',*  and  for  values  of  p  from  0.1  to 
0,5.  As  will  be  discussed  later,  the  empirical  limit  for  random  packing  of  monodisperse  spheres  is 
about  0.63,  but  in  practice  one  rarely  goes  much  beyond  about  0.5.  In  Figure  17,  we  note  that  the 
negative  teal  polarizability  manifests  itself  in  a  substantial  lowering  (diamagnetic  tendency)  of  p' 
from  unity.  In  Figure  18,  the  family  of  curves  for  p"  are  straightforward  reflections  of  the 
behavior  of  a"  (as  given  in  the  formula).  These  formulas  and  graphs  are  in  a  universal 
dimensionless  form,  expressed  in  tetmi  of  a  /  6.  It  remains  to  demonstrate  that  the  changing 
properties  are  relevant  to  the  materials  and  electromagnetic  parameters  of  this  investigation. 


IV.1.3  From  Model  Composite  to  Practical  Application 

At  this  point  we  transform  the  model  calculations,  expressed  in  dimensionless  ratios  of 
material  parameters,  into  teal  units  that  are  measurable.  The  major  parameter  is  the  resistivity, 
followed  by  particle  sizes  and  frequencies  of  interest  The  resistivity,  pf(,  for  the  alloy  Ni92Cr8  is 
55  x  10‘®  Am  [27].  When  this  is  substituted  into  the  formula  for  skin  depth,  6,  we  obtain  at  10 
GHz,  S  -  3.7  mm.  Proceeding  a  step  further  with  the  consideradon  that  peak  a"  (and  p")  occurs 
at  at  6*  2.3,  we  find  that  “important"  sphere  diameters,  dc  are  as  follows: 

</c-  17  pm  at  10 GHz 
de  “7.6  pm  at  50  GHz 
dc  -  5.4  pm  at  100  GHz 

Obviously  there  is  relevance  for  a  powder  with  particles  below  37  p  m  (-400  mesh). 

We  show  in  Figures  19  and  20,  the  calculated  model  behavior  for  p'  and  n"  for  this  alloy 
composidon,  for  a  volume  loading,  p  *  0.4,  and  for  various  sphere  diameters  of  interest,  i.e.,  5, 
10, 20,  and  30  p  m  (Similar  calculations  for  other  values  are  presented  below).  It  should  be 
noted  that  our  calculations  above  have  been  in  the  cgs  system  for  which  p\  p"  are  dimensionless. 
As  we  move  closer  to  real  electrical  measurements  in  the  “more  practical"  mks  SI  system,  these 
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Figure  17.  Magnetic  polarliablllty  (rail  part)  calculated  for  compoaltea  of  conducting 
•plurai  at  varloua  voluma  loadlnga,  p.  m'  va.  (a  16) 
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Figure  18.  Magnetic  permeability  (Imaginary  part)  calculated  for  compoaltee  of  conducting 
•pherea  at  varloua  volume  loadlnga,  p.  t u£  va.  (a/S)  . 
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Flgurt  19.  Relative  real  magnetic  permeability,  va.  frequency  calculated  for  compoiltee 
of  conducting  spberei  of  varloue  dlameterii  composite  loading  volume,  p  ■  0.4| 
sphere  alloy  resistivity,  pu-  55  x  10**  flm. 
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Figure  20.  Relative  Imaginary  magnetic  permeability,  vs.  frequency  calculated  for 
composites  described  In  Figure  19. 
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expressions  or  graphs  of  n\  \i*  (and  later  of  e\eH )  are  to  be  considered  as  relative  to  the 
permeability  (or  permittivity)  of  free  space.  In  other  words  we  have  n'r*i m'  /  ju,.  etc. 

The  qualitative  remarks  made  for  Figures  17  and  18  remain  true  for  these  figures. 
However,  the  separation  of  the  features  for  the  various  sphere  diameters  is  rather  dramatic.  It 
certainly  invites  experimental  demonstration  in  that  the  range  of  sizes  and  of  frequencies  for  this 
investigation  spans  the  calculated  behavior  very  nicely.  It  is  for  these  reasons  that  we  prepared 
samples  with  different  particle  sizes  (and  volume  loadings)  as  listed  in  Table  3;  i.e.,  -10  n  m, 

+  10-20  /im,  and  +  20-37  /urn, 

It  is  instructive  to  extend  the  critical  diameter/frequency  consideration  somewhat  further  on 
either  end.  If  we  move  to  the  far-infrared,  to  a  frequency  of  1  THz  (1  x  1012  Hz)  alternatively 
expressed  in  wave  numbers  as  33  cm'1,  the  critical  diameter  for  maximum  absorption  in  this  alloy 
would  be  1.7  fi  m,  which  is  in  line  with  predicdons  for  such  research  [25].  This  domain  is  of 
Interest  for  selective  absorbers  of  solar  or  infrared  radiation.  At  the  other  extreme  that  of  low 
frequency,  the  critical  diameter  at  10  kHz  is  1.7  cm.  This  is  the  regime  for  induction  heatingl 

In  each  of  these  two  limiting  regions  one  is  unlikely  to  find  results  obtained  and  presented 
in  terms  of  the  complex  permeability  (and  permittivity).  Instead  the  focus  is  on  absorption  per  se, 
because  of  measurement  techniques  as  well  as  because  of  practical  interest.  For  these  reasons, 
measurements  of  the  permeability  reported  here  make  a  useful  contribution  to  the  understanding  of 
artificial  dielectric  materials  in  the  microwave  frequency  region,  and  they  relate  to  the  nearby  far- 
infrared  region  also.  At  present,  we  are  aware  of  only  one  report  which  presents  these  kinds  of 
measurements  and  analysis.  It  deals  with  systems  of  millimeter-sized  ferromagnetic  steel  spheres 
over  the  kilohertz  and  megahertz  frequency  range  [28]. 

There  is  yet  another  interest  in  such  model  results,  a  technological  one.  The  RAM  physics 
model  [1]  predicts  that  high  frequency  loss  in  a  magneto-dielectric  should  become  negligible 
beyond  certain  frequencies  and/or  temperatures  (depending  on  material  parameters).  When  the 
induced  magnetic  effects  are  included,  one  may  expect  to  find  lossy  properties  outside  the  limits 
indicated  previously. 
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IV.2  Alternative  Approaches  To  Composite  Dielectrics 

The  model  dielectric  of  interest  herein  represents  one  limit  in  the  properties  of 
inhomogeneous  media  (or  composite  materials)  which  have  provided  problems  for  theorists  and 
engineers  for  over  150  years  [2].  It  is  relevant  to  note  that  for  composite  materials  the  calculations 
of  dielectric  constant,  magnetic  permeability,  electrical  conductivity,  thermal  conductivity  and 
diffusion  properties  can  often  be  formulated  in  precisely  analogous  ways.  Thus  there  is  an 
extensive  common  literature  for  such  phenomena  from  which  one  can  readily  apply  theoretical  and 
experimental  results  to  the  various  properties. 

There  is  a  basic  distinction  between  several  classes  of  models  of  heterogeneous  dielectric 
media.  More  modem  theoiy  follows  the  work  of  Bruggeman  [2, 29].  This  is  a  symmetrical  (two- 
component)  effective  medium  model  wherein  particles  of  type  1  are  mixed  randomly  in  some 
proportion  with  those  of  type  2,  so  that  each  may  be  surrounded  by  particles  of  both  types.  If  one 
type  is  conducting  and  the  other  insulating,  there  will  be  a  percolation  threshold  composition  at 
which  conducting  paths  are  established.  The  volume  loading  for  this  to  occur  is  about  0.15  to  0,2. 
This  model  is  often  labeled  EMA  (eff.  med.  approx.).  The  upper  inset  in  Figure  21  depicts  this 
topology,  and  the  curve  shows  the  dramatic  change  that  occurs  at  percolation  according  to  one 
mathematical  model. 

In  another  (older)  class  of  model  one  assumes  that  particles  of  type  1  are  completely 
surrounded  by  material  of  type  2.  It  is  clearly  necessary  to  specify  whether  metal  has  been  added 
to  a  dielectric  host  or  vice  versa,  i.e.,  the  model  is  unsymmetrical.  This  class  of  model  is 
associated  with  the  names  of  Maxwell,  Maxwell-Gamett,  Clausius  and  Mossotti,  and  Lorentz  and 
Lorenz.  It  too  is  really  an  effective  medium  (or  molecular  field)  model.  It  is  variously  called  by 
any  of  the  names  just  cited,  e.g.,  MGT,  CM,  etc.  Real-life  examples  of  this  model  are  ordered 
(cubical)  arrays  of  spheres  in  a  different  medium,  or  an  otherwise  correlated  medium  such  as 
coated  spheres  in  a  random  arrangement,  with  each  metal  particle  surrounded  by  a  concentric 
dielectric  coating.  The  lower  right  inset  in  Figure  21  describes  this  model,  and  the  curve  for  its 
behavior  is  indicated  Our  coated  alloy  particles  fall  into  this  claso  of  composites. 

In  recent  years,  rigorous  calculations  (and  some  measurements)  of  the  ordered  arrays  have 
appeared  [30, 31, 32, 33],  and  these  have  relevance  for  the  random  coaled  array.  Features  of  these 
theories  are  the  significant  effects  of  higher  order  multipoles  beyond  the  dipole-dipole  model  and 
the  (obvious  but  sometimes  ignored)  maximum  packing  limits,  Re-examination  of  the  earlier  data 
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Figure  21. 


Comparison  of  classical  dielectric  models  of  percolating  (Bruggeman)  vs.  non. 
percolating  (Maxwell/Clausius«Mossottl)  systems. 
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[1]  confirms  the  potential  applicability  of  these  recent  calculations,  and  deviations  therefrom  cast 
some  light  on  the  microstructure  of  the  random  coated  arrays. 


IV.2.1  The  Maxwell  (MGT,  CM)  Model 


An  excellent  derivation  of  the  formulas  for  this  model  with  emphasis  on  the  concepts  and 
approximations  is  presented  by  Landauer  [2].  A  contrasting  experimentally  oriented  discussion  on 
methods  of  approximating  the  behavior  of  two-phase  (heterogeneous)  systems  has  been  developed 
by  Mitoff  [21]. 

Mitoff  s  report  also  calls  attention  to  a  very  simple  method  described  by  Lichtenecker  [34] 
for  approximating  the  properties  of  two-phase  systems,  particularly  when  the  properties  of  the  two 
phases  are  not  very  different.  Although  the  method  is  not  fundamentally  correct,  it  is  easy  to  use 
and  does  agree  with  experimental  results. 

Lichtenecker's  rule  states  that  a  plot  of  log  k  (k  being  the  property  value  or  “constant”  for 
the  composite)  versus  vr '  me  fraction  p  of  phase  2  in  phase  1  is  a  straight  line  between  the  end 
values,  A  graphical  comparison  with  a  correct  treatment,  such  as  Maxwell's  (discussed  below) 
shows  that  the  difference  between  the  two  is  insignificant  with  respect  to  experimental  accuracy 
where  the  constants  of  the  two  phases  are  within  an  order  of  magnitude  of  each  other. 


Maxwell  derived  an  equation  which,  at  infinite  dilution,  exactly  describes  the  k  of  a  mixture 
of  spheres  with  property  constants  kj  and  k2,  ( k2>k] ).  While  he  used  kc  and  k<i  for  the  continuous 
and  the  discontinuous  phases,  Mitoff  recasts  this  into  two  expressions: 


k-k^ 


-p)(k,-kS 

*, +:*,+<!- f  )(*,-*,) 


for  ph°se  2  being  the  continuous  phase;  and 


ki+2ki-2p{ki~ki) 
Jhj  +  2  fcj  +  p(ki  “fcj) 
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for  phase  1  continuous.  It  assumes  that  the  discontinuous  phase  consists  of  spheres  and  that  the 
spheres  are  separated  enough  not  to  disturb  the  lines  of  “flux"  around  a  neighbor.  It  has  been 
shown  by  experiment  to  be  quite  accurate  up  to  10%  dilution  and  is  a  fair  approximation  in  some 
cases  even  up  to  50%. 


An  alternative  writing  of  the  second  of  the  above  equations  may  be  more  familiar  to  some. 
Thus  for  the  Maxwell,  Maxwell-Gamett,  Clausius-Mossotti  model  we  also  have 


k*k j 


i  »*(*-*> 

(1  -  p)*  +  (2  +  p)fej 


When  k2»kj  and p  is  small  (low-concentration)  this  formula  reduces  toward 

*■*(1+3  P) 


Although  equations  of  the  general  type  above  have  been  available  for  many  years,  they 
have  seen  little  use,  probably  because  they  are  not  the  kind  one  keeps  in  one's  head  for  a  quick 
approximation.  Mitoff  s  report,  however,  shows  that  a  very  simple  treatment  applies  when  the 
two  (separate)  pr  .ties  are  very  different.  This  is  precisely  the  region  for  which  the  simple 
Lichtenecker  method  does  not  apply, 

For  large  differences,  extending  to  a  number  of  orders  of  magnitude,  one  adopts  the 
procedure  of  plotting  log  k/kj  versus  volume  fraction.  From  Mitoff,  we  obtain  Figure  22, 

presenting  the  Maxwell  curves  for  kilki  ratios  from  100,3  to  106,  The  upper  curves  represent 
phase  2  continuous  and  the  lower  ones  phase  1  continuous.  The  straight  line,  L,  represents  the 
Lichtenecker  approximation.  One  notices  that  when  kitki  is  less  than  10,  the  mathemadcally 
defensible  Maxwell  treatment  is  not  much  different  from  Lichtenecker's. 


Mitoff  goes  a  step  further  in  observing  that  the  limiting  behaviors  in  Figure  22  for  large 
differences  (&2»&;)  are  really  straight  lines  with  slopes  that  do  not  involve  the  property  constant 
of  the  discontinuous  phase.  He  obtains  the  relationships: 


Alog10(*/*) 
A(1  -P) 


a  -0.65  (*  ■ 


and 
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1.3  (*,«*,) 

Ap 


These  are  shown  graphically  on  Figure  23.  So  long  as  the  It's  are  very  different,  one  expects  the 
points  to  fall  on  these  straight  lines,  especially  at  volume  fractions  close  to  the  kc  side. 

Mitoff  compares  this  approximate  method  with  the  complete  Maxwell  equation  and  with 
Lichtenecker's  rule  for  a  wide  range  of  ratios  of  k2lkj,  He  shows  that  for  one-order-of-magnitude 
difference  or  more  the  first  or  upper  straight  line  for  ki  »  kc  is  useful  even  out  to  rather  high 
volume  fractions  of  the  discontinuous  phase.  At  less  than  an  order-of-magnitude  difference, 
Lichtenecker's  line  is  a  good  approximation  to  the  Maxwell  equation  while  the  others  are  not. 

Inspection  of  Figure  22  shows  that  the  Maxwell  (MOT,  CM)  expressions  continue 
analytically  across  the  whole  range  of  volume  fraction,  p.  That  the  curves  are  shown  dashed  at 
extreme  p  values  is  a  recognition  of  two  points:  one  is  the  assumption  of  non-contact,  translated  as 
meaning  a  dilute  system;  the  second  is  the  fact  that  the  volume  loading  of  uniformly  sized  spheres 
in  random  packing  cannot  exceed  0.63,  and  reaches  this  only  under  specially  favorable  conditions. 
We  return  to  this  point  below.  Nevertheless  the  Maxwell  (MOT,  CM)  formula  is  useful,  as  is  the 
Mitoff  approximation  when  used  with  care. 

There  are  two  applications  of  practical  interest  for  this  discussion.  The  upper  branch 
curves  of  Figure  22  (or  their  analytic  expression)  successfully  describe  the  effects  of  porosity  in  a 
host  medium.  The  second  application  of  special  interest  is  that  of  a  coated  sphere.  It  connects  to  a 
designated  task  of  Table  1,  i.e.,  “Analyze  the  effects  of  microstructure  control  in  passivated  metal- 
insulator  composites." 

In  Section  Q.5  we  examined  the  surface  condition  of  the  metal  powder  and  concluded  that 
there  is  an  insulating  oxide  coadng  of  about  0.1  to  0.2  p  m  thickness.  It  would  be  rash  to  call  this 
a  passivating  layer,  however.  This  conclusion  derives  from  the  considerable  reaction  suffered  by 
this  alloy  powder  when  composite  preparation  using  an  inorganic  binder  was  attempted,  as 
discussed  in  Section  in. 

The  presence  of  the  insulating  coating  is  very  significant  for  the  electromagnetic  behavior, 
especially  for  the  dielectric  properties  (See  Figure  21).  Loosely  speaking,  it  prevents  the  individual 
metal  particles  from  shorting  to  each  other,  More  elegantly,  it  means  we  have  a  correlated 
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Figure  22.  Maxwell  (MOT,  CM)  equations  for  mixtures,  on  a  semilog  plot.  Upper  set  ti  for 
Phau  2  continuous |  lower  set  Is  for  Phase  1  continuous;  L  1s  Llchtenecker's 
approximation.  (The  number  on  each  curve  signifies  the  power  of  ten  In  the 
ration  k2/ki.)  (From  Reference  21) 
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Figure  23.  MltofPs  graphical  approximations  to  the  Maxwell  (MGT,  CM)  equations  for 
k2  »  ki.  (From  Reference  21) 
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heterogeneous  medium  in  which  conducting  particles  are  completely  surrounded  by  insulating 
material.  This  restricts  the  theoretical  approach  to  describing  the  dielectric  behavior.  Volume 
loadings  up  to  or  approaching  the  packing  limits  can  be  achieved  with  coated  particles  (spheres) 
without  the  establishment  of  conducting  paths,  i.e.,  percolation. 

With  a  coating  on  the  conducting  sphere,  as  well  as  a  binder  material,  the  resulting 
composite  is  a  three-phase  medium.  However,  it  is  possible  to  replace  the  coated  sphere 
combination  by  a  single,  slightly  modified,  dielectric  sphere,  at  least  in  our  frequency  regime.  An 
apparently  direct  application  of  the  Maxwell  (MOT,  CM)  expressions  yields  [35]  the  single 
dielectric  constant 


2c,  +  e1  +  2fl3(gi  -£,)' 

2£,  +  £J-Q,(C2-Ct) 


where 


with  t  the  coating  thickness  and  a  the  sphere  radius.  For  £,  »  c, ,  medium  1  being  the  insulating 
coat,  this  reduces  to 


«£, 


‘1+2  Qr 

„1-GS. 


For  our  estimated  coating  thickness  of  0.1  n  m,  on  a  sphere  of  average  diameter  23  /.t  m,  we  find 
£«/£,  «  115.  For  simplicity  we  take  £,  for  the  coating  as  that  for  Cr203,  i.e.,  11.1,  from  which 
ea  ■  1276.  This  value  for  £„  is  still  so  much  greater  than  that  for  typical  binders  that  the  model 
calculations  are  unaffected 


IV.2,2  The  Bruggeman  Effective  Medium  Approximation 

We  described  the  initial  elements  of  this  EMA  model  in  IV.2  above.  A  visual  image  is 
presented  in  Figure  21  (after  Landauer  [2]).  The  two  components  are  mixed  randomly.  The 
shaded  crystallite  of  type  1  (of  volume  fraction,  p)  is  surrounded  by  crystallites  of  both  types.  In 


59 


the  analysis  [2],  the  surrounding  (two-phase)  material  is  imagined  as  a  single  medium  of  uniform 
properties.  The  homogeneous  dielectric  constant  eg  in  this  case  arises  from  the  solution  of  a 
quadratic  equation  which,  assuming  spherical  particles,  is: 


2ea  +  £, 


+(i-p)tiz£2—o 


2ea  +  e2 


At  quite  low  co.  Titrations,  p,  the  EMA  and  the  Maxwell-CM  give  identical  results.  The  EMA  has 
the  appealing  ar  Jt  of  treating  panicles  of  types  1  and  2  on  an  equivalent  basis  (hence  a 
symmetrical  theory). 

As  we  have  discussed  above,  our  artificial  dielectric  is  one  with  a  correlation  between  the 
two  phases  of  the  medium.  This  rules  against  the  use  of  the  EMA  model  and  in  favor  of  one  for 
which  the  metal  is  embedded  in  the  insulating  medium  (thanks  to  the  coating). 

An  interesting  and  distinctive  feature  of  the  EMA  is  that  it  predicts  a  metal-insulator  or 
percolation  transition.  In  this  version  it  occurs  atp»  1/3,  while  in  other  theories  this  value  can  be 
as  low  as  p»  0.15.  (For  quasi-two-dimensional  mixtures,  the  values  are  somewhat  higher).  A 
good  illustration  of  this  prediction  experimentally  is  shown  in  Figure  24  [36].  For  these  samples 
the  critical  percolation  fracdon  (vol.  loading)  is  about  0.20.  This  demonstration  is  a  convincing 
argument  against  the  EMA  class  of  model  for  the  artificial  dielectrics  of  this  study. 


IV.2.3  Exact  Calculations  for  Arrays  of  Spheres 

Concurrent  with  the  interest  in  random  heterogeneous  media  is  an  attention  to  the  equivalent 
problem  for  ordered  arrays  of  (conducdng)  elements.  Lord  Rayleigh  in  1892  showed  a  significant 
advance  beyond  Maxwell's  work  when  he  took  into  account  induced  octupole  moments.  The 
ordered  array  lends  itself  to  a  more  rigorous  formulation  of  the  problem.  The  subject  received  a 
technical  boost  with  the  suggestion  of  Kock  [37]  that  artificial  dielectrics  (of  the  ordered  type) 
could  be  used  as  microwave  lenses.  There  has  since  been  a  great  deal  of  engineering  interest, 
involving  a  variety  of  conductor  shapes  [38].  We  have  found  particularly  useful  a  paper  by 
Kharadly  and  Jackson  [30]  both  for  the  (lossless)  effective  permittivity  calculations  and 
measurements,  and  for  discussions  of  the  dielectric  loss.  Most  recently,  with  modem  computers, 
one  can  carry  out  numerical  calculations  of  these  problems  to  arbitrary  precision  limits  which  were 
totally  unattainable  heretofore  [31-33, 39,40]. 
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Figure  24.  The  dielectric  constant  of  several  Ag-KCI  samples  aa  a  function  of  metal  volume 
fraction,  p.  (From  Reference  34) 


Figure  2S.  Dielectric  constants  for  cubic  lattices  as  calculated  In  Reference  31  (upper)  and 
References  32,  33  (lower) 


As  noted  above,  the  higher  order  multipoles  (beyond  the  dipole  terms  of  the  Maxwell/CM 
theories)  have  an  important  effect  as  the  spherical  conductors  get  closer  at  increased p.  A  measure 
of  the  success  of  the  calculations  is  that  the  calculated  permittivities  exhibit  an  anomaly  at  the 
critical  volume  loading  where  the  spheres  just  touch,  in  contrast  with  the  Maxwell  (MOT,  CM) 
models.  The  conceptual  drama  of  this  feature  shows  in  the  simplicity  of  a  confirmatory  experiment 
for  the  body  centered  cubic  lattice  [33],  It  suffices  to  construct  only  a  one-eighth  portion  of  a  unit 
cell  using  conducting  and  insulating  plates  plus  two  one-eighth  sections  of  metal  ball  bearings! 

The  metal-insulator  "anomaly"  or  transition  when  they  touch  is  thus  lattice  dependent.  This 
follows  from  the  3-dimensional  packing  limits  for  hard,  uniform  spheres,  as  follows, 
simple  cubic  p  ■  it  /  6  ■  0.32 
body  centered  cubic  p  «  rrV3  /  8  «  0. 68 
face  centered  cubic  p  *  «V2  /  6  *  0.74 
or  hexagonal  close-packed 

In  Figure  23  we  reproduce  the  results  of  the  calculations,  the  upper  from  Doyle  [3 1,39]  and 
the  lower  from  the  Australian  group  [32,33,40].  The  Maxwell-CM  curve  in  the  upper  portion 
agrees  with  the  exact  calculations  out  to  p  m  0.33  for  SC  and  to  p  a  0.3  for  BCC  and  FCC.  This 
provides  the  justification  for  the  reasonable  success  in  our  prior  work.  It  was  the  right  type  of 
model  and  a  not-bad  approximation.  On  the  other  hand,  the  Maxwell/CM  model  fails  to  diverge 
until  p  tends  toward  unity.  The  line  marked  R  is  the  Rayleigh  calculation  for  the  SC  lattice 
(corrected  [31])  showing  some  multipole  effects  but  failing  to  diverge  at  the  contact  limit.  In 
addition  experimental  data  are  included  for  the  SC  lattice  which  are  in  line  with  the  exact  results. 

The  calculations  of  the  Australian  group  in  the  lower  portion  of  Figure  25  agree  in  detail 
with  those  of  Doyle,  since  the  methods  of  attack  were  equivalent.  They  have  added  another  curve 
estimated  by  scaling  arguments  and  intended  to  apply  to  a  random  or  disordered  array.  The 
empirical  packing  limit  in  this  case  (as  noted  above)  is  p  ■  0.63,  although  that  limit  is  tricky  to 
reach  and  values  nearer  top  *  0.60  are  more  likely.  They  even  express  this  curve  with  a  closed 
approximate  formula  as  follows: 

e  -  0.05398  -  2.0673  x  ln(0. 63  -  p). 
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IV, 2.4  Reanalysis  of  Prior  e'  Results 

In  our  prior  study  [1]  we  obtained  dielectric  measurement  data  on  a  composite  system  of 
magnetic  metal  particles  in  the  same  polyurethane  binder  as  employed  herein  (we  call  this  the  B- 
series).  The  particles  were  reasonably  spherical  and  had  an  insulating  coating  as  manifested  by  a 
resistance  test  (cf.  the  last  entry  in  Table  2).  The  composites  studied  had  (magnetically  measured) 
volume  loading  ranging  from  about  p  *  0.2S  to  p  ■  0.50.  It  is  instructive  to  re-examine  those  data 
in  the  light  of  the  preceding  discussion,  especially  as  regards  exact  calculations  for  arrays, 
inasmuch  as  the  prior  analysis  had  considered  only  the  Maxwell  (MOT-CM)  model. 

In  Figure  26  we  display  that  e *  data  in  comparison  with  the  Maxwell-Mitoff 
approximation.  The  various  data  symbols  refer  to  particle  groups  with  differing  average  sizes  and 
size  ranges.  Most,  but  not  all,  of  the  groups  have  data  that  are  approximated  by  the  solid  line 
drawn  as  a  guide  to  the  eye.  (We  cannot  really  explain  that  “maverick"  group,  but  we  note  that  it 
was  comprised  of  the  widest  range  of  particle  sizes  of  the  set,  about  a  factor  of  four  from  the 
largest  diameter  to  smallest.)  We  confess  that  in  the  prior  analysis  [1],  we  were  so  infatuated  by 
the  potential  of  the  Mitoff  approximation  line  (of  predetermined  slope-  see  Section  IV.2.1)  that  we 
used  that  line  fitted  to  the  lower  p- values  of  the  B-serics  data  set.  (Those  p-values  are  not  really 
very  low!)  In  that  case,  in  the  limit  of  p  ■  o,  the  line  came  close  to  but  still  deviated  ftom  the 
separately  measured  binder  value,  e, .  In  the  greater  wisdom  of  hindsight,  Figure  26  fixes  the 
Mitoff  line  to  originate  at  the  binder  value.  We  strongly  infer  ftom  this  presentation  that  upward 
deviations  ftom  the  Maxwell/Gausius-Mossotti  model  occur  experimentally  at  lower  values  of 
volume  loading,  p,  than  we  suspected  earlier. 

Finally,  we  offer  a  comparison  of  the  experimental  representation  with  some  of  the  better 
theoretical  models.  This  is  shown  in  Figure  27  where  the  dashed  line  portrays  the  behavior  of 
most  of  the  B-series,  apart  ftom  the  maverick  set.  All  are  plotted  as  a  ratio  to  the  binder  e't  value. 

One  comparison  is  to  the  predicted  behavior  for  a  simple  cubic  array  (S.C.)  calculated  on  a 
rigorous  model  [31, 32],  which  is  the  array  diverging  at  the  lowest  p-value  (p  -  0.523).  The 
second  comparison  is  to  the  quasi-empirical,  quasi-theoretical  model  for  a  nonpercolating  random 
composite  [33].  Ibis  idealized  random  model  was  constructed  to  diverge  at  the  empirical  packing 
limit  for  random  systems,  p  si  0.63  and  fit  to  a  logarithmic  divergence  that  appears  to  describe  the 
rigorously  calculated  array  models. 
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The  data  curve  mimics  the  upward  curvature  of  the  better  theoretical  models  but  precedes 
their  predictions  noticeably,  In  discussions  with  our  consultant,  Professor  W.T,  Doyle  of 
Dartmouth  College,  he  suggested  an  attractive  explanation  based  on  the  marked  non-linearity  of  the 
€ '  vs.  p  curves.  In  contrast  to  the  lattice  arrays  of  particles,  one  could  have  fluctuations  of  volume 
loading  in  a  randomly  packed  composite.  Indeed  such  is  the  norm  in  random  systems,  but  we 
imagine  that  on  some  spatial  scale  it  can  make  a  difference.  (Note:  Our  use  of  semilogarithmic 
coordinate  paper  hides  some  of  the  non-linearity.)  It  is  reasonably  clear  that  there  is  room  for 
improved  modeling. 


IV, 3  Improved  Modeling  of  Randomly  Packed  Artificial  Dielectrics;  Pair 

Interaction  Model 

In  discussion  with  W.T.  Doyle,  we  conceived  of  an  instructive  but  very  much  simplified 
approach  to  the  problem  of  randomly  packed  artificial  dielectrics.  The  germ  of  the  idea  comes  from 
inspection  of  a  micrograph  of  one  of  our  actual  composites  of  Ni-Cr  particles  in  polymer  at  10% 
volume  loading  (Figure  10a),  in  comparison  with  a  schematic  diagram  used  by  mlcroscoplst* 
ceramists  to  estimate  porosity  shown  in  Figure  28.  In  a  limited,  but  incorrect  sense,  the  schematic 
diagram  may  be  termed  an  "Idealized  random"  composite  with  the  volume  loading  as  indicated 
While  there  is  no  periodicity,  it  lacks  the  density  fluctuations  of  a  truly  random  loading.  On  the 
positive  or  idealized  side,  the  particles  (■  dots)  never  touch  or  coalesce  until  the  packing  limit  is 
reached!  This  would  be  about  63%,  for  a  3-D  random  array.  Thus  this  "random"  schematic 
actually  represents  a  correlated  medium.  In  fact,  such  an  assumption  of  a  correlated  medium  was 
included  [41]  in  the  "estimated  calculation"  for  a  random  or  disordered  array  carried  out  by  the 
Australian  group  [33]  by  using  scaling  arguments  on  the  rigorous  results  for  an  ordered  array. 

In  contrast,  when  one  examines  the  micrograph  of  the  actual  composite  (Figure  10a),  local 
clustering  is  quite  evident.  This  comparison  led  us  to  consider  a  model  based  on  pair  interactions. 
Professor  Doyle's  report  on  this  model,  performed  as  a  subcontract,  follows. 


IV.3,1  Introduction 

The  average  permittivity  of  an  inhomogeneous  medium  composed  of  spheres  of  one 
substance  embedded  in  a  continuum  of  a  different  material  (cermet  topology)  is,  in  general,  a 
complicated  function  of  the  permittivities  of  the  constituents,  the  volume  filling  factor,  particle 
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Figure  21.  Schematic  diagram  of  "Idealised  random"  compoiltaa  adapted  rrom  figure  uied  by 
mlcroacopleta  to  eatlmate  poroalty  (For  comparlaoa  with  actual  composite  In 
Figure  10a.) 
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shape  and  size,  and  the  details  of  the  particle  distribution  function.  Although  the  permittivities  may 
be  complex  and  depend  upon  frequency,  we  shall  be  concerned  only  with  perfectly  conducting 
particles  in  the  electrostatic  limit 

In  the  simplest  case  of  a  regular  array  of  spheres  of  uniform  size  the  permittivity  may  be 
computed  exactly  [31-33].  When  one  or  the  other  of  the  various  parameters  has  a  non-zero 
variance,  the  problem  becomes  substantially  more  difficult  If  all  parameters  exhibit  wide 
variances,  the  problem  is  essentially  intractable.  Fortunately,  in  practice,  it  is  possible  to  obtain 
nearly  mono-disperse  distributions  of  good  spheres.  However,  even  an  array  of  uniform  sized 
spheres  presents  difficulties  when  the  spheres  are  randomly  distributed  in  space.  For  such  media 
the  permittivity  may  be  calculated  only  approximately. 

If  only  dipole  interactions  are  present,  the  permittivity  of  an  isotropic  medium  is  given  by 
the  well-known  Clausius  -Mossotti  relation.  For  regular  arrays  this  case  occurs  in  the  limit  of  low 
filling  factors.  For  random  distributions  close  encounters  can  occur  even  with  low  filling  factors, 
and  corrections  may  be  necessary  even  with  low  particle  densities.  Experiments  with  nominally 
“random”  media  often  exhibit  substantial  departures  from  the  predictions  of  the  approximate 
formulas,  usually  of  the  kinds  to  be  expected  from  enhanced  particle  interactions  due  to  particle 
clustering. 


IV.3.2  Clustering  Effects 

In  a  neutral  dielectric  medium  only  dipole  and  higher  multipolar  Interactions  exist,  and  all 
depend  inversely  upon  increasing  powers  of  the  distance  between  the  particles  of  the  medium.  As 
a  consequence,  the  multipole  interactions  grow  rapidly  with  increasing  volume  filling  factor. 
Because  they  are  positive,  all  multipole  contributions  lead  to  an  increase  in  the  dielectric  constant 
with  decreasing  interparticle  distance.  A  macroscopic  manifestation  of  this  is  a  strong  non-linearity 
of  the  dielectric  constant  when  plotted  as  a  function  of  volume  filling  factor.  All  experiments  and 
physically  valid  theoretical  approximations  show  this  nonlinear  behavior.  An  immediate 
consequence  of  this  simple  arithmetic  fact  is  that  any  departure  from  strict  constancy  of  the  density 
of  the  medium  leads  to  an  enhancement  of  the  dielectric  constant. 

Density  fluctuations  are  inescapable  in  a  random  medium.  They  occur  statistically,  if  for  no 
other  reason,  although  there  are  often  stronger  variances  associated  with  fabrication,  such  as 
incomplete  stirring,  panicle  cohesion,  settling,  substrate  distortion,  interparticie  forces,  and  other 
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factors.  Any  or  all  of  these  density  fluctuations  cart  be  expected  in  any  given  sample,  and  all  can 
lead  to  an  increase  in  the  measured  macroscopic  average  dielectric  constant  We  shall  refer  to  all 
such  enhancements  as  “cluster  effects." 

The  problems  associated  with  an  unrestricted  density  variation  are  formidable,  if  not 
actually  unsolvable.  In  this  report  we  shall  examine  one  of  the  simplest  types  of  clustering  effect: 
the  permittivity  enhancement  produced  by  the  presence  of  a  population  of  close  pairs  of  particles  in 
a  medium  consisting  of  perfectly  conducting  monodisperse  spheres.  We  shall  only  consider 
interactions  between  isolated  pain  of  particles.  All  larger  clusten  of  panicles  and  large-scale 
density  fluctuations  will  be  ignored, 

The  interaction  of  an  isolated  pair  will  be  treated  exactly,  including  the  effect  of  all  higher 
multipoles  mutually  induced  by  the  partners  on  each  other.  As  will  be  seen,  the  pair  Interaction 
leads  to  on  enhancement  of  the  effective  polarizability  of  the  pair.  Once  the  pair  polarizability  has 
been  found,  the  pair  may  be  treated  as  a  molecular  dipole  in  calculating  the  dielectric  constant  of  an 
ensemble  of  mixed  isolated  particles  and  pairs. 


IV.3.3  Pair  Polarizabilities 

Laplace's  equation  is  separable  in  bispherical  coordinates,  so  the  problem  of  two 
conducting  spheres  in  a  uniform  external  field  may  be  solved  exactly.  Using  this  coordinate 
system  Levine  and  McQuarrie  [42]  obtained  series  expressions  for  the  total  polarizability  of  a 
system  of  two  spheres  with  arbitrary  separation.  We  follow  their  treatment  here. 

The  polarizability  of  the  pair  may  be  written 

OW-  a,  +  «,  +  «„  (1 

where  a,  *  <xs »  R 1  is  the  polarizability  of  either  of  the  two  identical  isolated  spheres,  and  al2  is 
the  incremental  polarizability  caused  by  each  sphere's  udfil  field  on  the  other.  Since  die  solution  is 
exact,  this  enhanced  polarization  includes  the  effect  of  fill  higher  multipoles  induced  on  these  two 
spheres  in  a  uniform  external  field. 
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The  polarizability  tensor  of  a  pair  is  axially  symmetric,  so  there  are  two  independent 
components  a12tL  and  a12T(  longitudinal  and  transverse  to  the  pair  axis,  respectively.  The  average 
polarizability  of  an  ensemble  of  pairs  with  axes  distributed  isotropically  is  given  by 

**12, o  =  (**12  j,  ■*"  2i7-)  *  3  (2) 

For  an  isotropic  random  medium  a12i0  is  the  important  quantity. 

The  polarizabilities  a12  r  and  ali  L  are  given  by 


and 


«u,r  ■  2/?*sinh*  ^(-lr^sinhsTj,)"3 
#■  2 


(3) 


an,L  *4/?3sinh3  (4) 

where  rj0  is  defined  by 


(5) 


The  pair  separation  parameter  x  ■  r  /  2/?  is  determined  by  the  intetparticle  center-to-center 
distance,  r,  and  the  sphere  radius,  /?.  The  functions  are  given  S,(tj,),S, (t|#),  and  52(tj„)  by 


S 


A  (2r  +  l)* 

(To  exp[(2f  + 1)7?,]  - 1 


(6) 


The  polarizabilides  obtained  with  Equations  (2)  through  (4)  (normalized  to  ft’)  are 
tabulated  in  Table  8  and  shown  graphically  in  Figure  29.  From  the  figure  it  is  clear  that  a 
substantia)  enhancement  of  the  dielectric  constant  can  occur  in  the  presence  of  many  closely  spaced 
pairs.  For  spheres  Li  contact,  all  of  the  incremental  polarizabilities  are  of  the  same  order  of 
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magnitude  as  the  polarizability  of  an  isolated  sphere.  The  average  polarizability,  a12|0,  which 
determines  the  dielectric  constant  of  a  random  distribution,  drops  off  very  rapidly  with  increasing 
separation  of  the  spheres  in  a  pair. 


IV.3.4  Dielectric  Constant 

The  values  given  in  Table  8  art  for  pairs  of  particles.  The  incremental  polarizabilities  per 
particle  art  half  as  large  as  the  tabulated  values.  To  find  the  average  effective  polarizability  of  a 
particle  of  the  medium  we  must  add  one  half  of  the  incremental  polarizability  for  a  given  pair, 
ai2,o  U)/2' with  separation  parameter  x<  to  the  isolated  sphere  polarizability,  R\  and  then  average 
over  all  pairs  in  the  medium.  Thus 

(a)  ■  J?1  + 1  /  2j  «a(^)«i2,o4^  (7) 

where  *  Nt  /  N  is  the  fraction  of  pardcles  that  are  members  of  pairs,  X  is  the  interparticle 
distance  parameter,  and  R  is  the  particle  radius. 

Equation  (7)  shows  that  we  need  Anther  information  about  the  particle  distribution  in  order 
to  calculate  the  effective  polarizability.  Tills  kind  of  detailed  structural  knowledge  is  required  with 
any  approximation  higher  than  the  dipole  approximation.  In  the  present  case  we  must  know,  or 
make  some  assumption  about,  the  number  of  pairs  present  for  each  value  of  the  separation 
parameter  £  ■  2 R, 

The  simplest  assumption  we  can  make  regarding  ^(x)  is  that  all  of  the  pairs  have  the  same 
separation,  (%),  say.  Then  we  have 

<o)  -  J?*[l+ nJCr)o,I.0(je)/2«*]  (8) 

so  that  there  are  only  two  adyustable  parameters:  the  fractional  number  of  particles  in  pairs,  ^(x), 
with  the  effective  separation  parameter  £.  The  value  of  aia,0(7) is  taken  from  Table  8,  using  the 
appropriate  average  particle  separation  parameter  %, 
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Table  8 

INCREMENTAL  POLARIZABILITIES,  al2,  PER  PAIR  (NORMALIZED  TO  R^) 
AS  A  FUNCTION  OF  THE  SEPARATION  PARAMETER,  x  ■  r/2R 

Subscript!  T,  L,  and  O  designate  the  transverse,  longitudinal 
'  and  average  composites,  respectively.  (Sphere  radius  Is  R,  and 
Interparticle  separation,  center-to-center,  is  r.) 


a 

®U,L 

«U,0 

5.00000 

-0.0020 

0.0040 

0.0000 

4.00000 

-0.0039 

0.007B 

0.0000 

3.00000 

-0.0092 

0.0187 

0,0001 

2.00000 

-0.0307 

0,0648 

0.0012 

1,90000 

-0.0354 

0.0762 

0.0016 

i.eoooo 

-0,0417 

0.0904 

0.0023 

1.70000 

-0.0492 

0.1086 

0.0034 

1.40000 

-0.0SS4 

0.1323 

0.0050 

1.50000 

-0.0704 

0.1641 

0.0078 

1.40000 

-0.0853 

0.2083 

0.0126 

1.30000 

-0,1043 

0.2729 

0.0215 

1.20000 

-0.1295 

0.3737 

0.0396 

1.10000 

-0.1591 

0.5695 

0.0838 

1.09000 

-0.1425 

0.5998 

0.0916 

1.00000 

-0.1440 

0.6334 

0.1005 

1.07000 

-0.1494 

0.6719 

0,110? 

1.04000 

-0.1733 

0.7159 

0.1231 

1.05000 

-0.1770 

0.7673 

0.1377 

1.04000 

-0.1809 

0.8291 

0.1558 

1.03000 

-0,1848 

0.9066 

0.1790 

1.02000 

-0.1887 

1.0107 

0.2111 

1.01000 

-0.1928 

1,1730 

0.2625 

1.00900 

-0.1932 

1,1959 

0.2698 

1.00800 

-0,1934 

1.2210 

0.2779 

1.00700 

-0.1940 

1,2487 

0.286? 

1.00400 

-0.1944 

1.2797 

0.2970 

1.00500 

-0.1948 

1.3152 

0.3085 

1 , 00400 

-0.1953 

1.3569 

0.3221 

1.00300 

-0,1957 

1.4078 

0.3388 

1.00200 

-0.1941 

1.4746 

0.3608 

1.00100 

-0.1945 

1,5766 

0.3945 

1.00090 

-0.1945 

1,5908 

0.3992 

l.oooeo 

-0.1944 

1.6064 

0.4044 

1.00070 

-0.1944 

1,6236 

0.4101 

1  .00040 

-0,1947 

1.6429 

0.4165 

1 .00050 

-0.1947 

X  *  6630 

0.4239 

1.00040 

-0.1947 

1.6910 

0.4325 

1.00030 

-0.1948 

1.7229 

0.4431 

1 .00020 

-0.1948 

1,7649 

0.4571 

1.00010 

-0.1949 

1,1299 

0.4787 

1 .00009 

-0.1949 

1.8391 

0,4818 

1 .00009 

-0.1949 

1.8492 

0,4851 

1.00007 

-0.1949 

1,8604 

0.4889 

1 .00004 

-0.1949 

1 .8730 

0.4931 

1.00005 

-0.1949 

1.8874 

0.4979 

1.00004 

-0.1969 

1.9046 

0,5036 

i*.  00003 

-0.1949 

1.9257 

0.5106 

1 .00002 

-0.1949 

1.9540 

0,5200 

1 .00001 

-0.1949 

1.9983 

0.534B 

1 .00000 

-0.1969 

2.8082 

0.8048 
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The  effective  polarizability  is  then  used  in  the  dipole  approximation  for  the  dielectric 
constant.  The  Clausius-Mossotti  expression 


n  '♦**<«> 

l-ijrW(a) 


is  consistent  With  our  assumption  that  only  pairs  and  isolated  spheres  are  present. 


(9) 


IV.3.5  Comparison  with  Prior  Experiment 

Our  results  apply  to  a  suspension  of  conducting  spheres  in  the  long  wavelength 
(electrostatic)  limit  with  intermediate  rilling  factors.  Although  the  Clausius-Mossotti  expression  for 
isolated  spheres  is  independent  of  the  sphere  size,  the  pair  calculations  assumed  equal-sized 
spheres,  so  we  should  restrict  our  comparison  to  fairly  mono-disperse  samples.  Some  of  the  older 
literature  reviewed  by  DcLoor  [43]  report  dielectric  constant  enhancements  that  can  be  attributed  to 
cluster  effects.  The  measurements  of  Guillien  [44]  on  emulsions  of  Hg  in  oil  are  of  particular 
interest,  because  both  particles  and  host  are  liquids,  so  that  the  particles  are  exceptionally  smooth 
and  spherical.  We  shall  first  compare  our  calculations  with  Guillien's  low  frequency  results  and 
with  the  microwave  measurements  reported  in  our  prior  study  [1]. 

The  largest  possible  isotropic  dielectric  constant  enhancement  obtainable  with  isolated  pairs 
occurs  when  all  of  the  spheres  form  pairs  in  contact  The  number  of  isolated  spheres  is  then  zero, 
and  the  number  of  pain  is  maximum,  provided  that  no  sphere  is  permitted  more  than  a  single  point 
of  contact.  Prom  Table  8  we  see  that  at  contact  the  average  polarizability  per  sphere  is  equal  to 
1.4024  R3.  The  solid  lines  in  Figures  30  and  31  show  the  Clausius-Mossotti  results,  with  and 
without  pairing.  The  bottom  curve  is  the  Clausius-Mossotti  result  with  no  pain  present.  The 
upper  two  curves  show  the  pair-enhanced  dielectric  constant  according  to  two  extreme  assumptions 
about  the  dependence  of  the  pairing  probability  upon  particle  density,  as  discussed  next. 

In  the  topmost  curve  of  Figures  30  and  31  we  assume  that  pairing  probability  is 
independent  of  particle  density,  i.e.,  n,  -  l.  One  conspicuous  consequence  of  this  assumption  is 
that  the  derivative  of  the  curve  of  dielectric  constant  versus  rilling  factor  is  larger  than  the  Clausius- 
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Figure  30. 
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VOLUME  LOADING,  p 

Dielectric  constant.  «r',  va.  volume  loading,  p,  Calculated  modala  with  and 
without  pairing,  ng.  Experimental  data  from  Reference  1. 
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VOLUME  LOADING,  p 

Figure  31.  Reduced  permittivity.  ve.  volume  loading,  p.  Calculation!  ai  for  Figure 

30.  Data  for  iphtrlcal  drops  of  Hg  In  heavy  oil,  from  Reference  44. 
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Mossotd  value  3et  at  p  ■  0.  This  is  admittedly  an  extreme  assumption,  in  that  complete  pairing  is 
assumed  even  in  the  limit  of  very  small  filling  factor.  Such  a  high  clustering  tendency  may  be 
appropriate  for  cohering  particles  that  resist  breakup  by  stirring.  Some  of  the  experimental  data 
seem  to  involve  some  low  density  clustering.  Whether  or  not  such  media  exist,  the  example  serves 
to  show  that  even  clusters  as  small  as  pairs  can  produce  significant  effects. 

Another  simple  case  is  obtained  by  assuming  that  the  pair  density  is  simply  proportional  to 
the  filling  factor  p.  Assuming  that  n^x)  ■  P>Pc  ■  we  have 

<«)-*’[i+o,ua)P/2p./e]  do) 

Setting  pc  -  0.63,  the  value  appropriate  for  a  random  medium,  we  get  the  middle  curves  in  Figures 

30  and  31.  This  is  a  more  natural  assumption  for  noncohering  particles,  since  for  such  materials 
one  would  expect  the  pair  density  to  go  to  zero  in  the  limit  of  low  filling  factors.  With  this 
assumption  there  are  no  pain  present  at  zero  filling  factor,  so  the  initial  slope  is  equal  to  3e,  in 
agreement  with  the  Clausius-Mossotti  result. 

The  calculated  curves  are  the  same  in  Figures  30  and  31  except  for  normalization.  The 
points  shown  are  the  experimental  results,  Figure  30  shows  microwave  measurements  on  solid 
particles  from  Reference  1.  Microscopic  examination  showed  the  samples  to  be  relatively  mono- 
disperse  distributions  of  nearly  spherical  particles.  Figure  31  shows  Ouillien's  measurements  on 
Hg  emulsions  in  oil  [44],  The  similarity  of  the  experimentally  measured  dielectric  enhancements  of 
the  solid  samples  in  Figure  30  and  the  liquid  sample  in  Figure  31  is  worth  noting,  since  the  Figure 

31  liquid  sample  presumably  exhibits  an  enhancement  produced  by  density  variation  alone, 
uninfluenced  by  a  distribution  of  particle  shape. 

IV.3.6  Conclusions  (Doyle  Subcontract) 

Using  the  polarizabilities  in  Table  8,  together  with  Equations  (7)  and  (9),  one  could  easily 
calculate  the  dielectric  constant  in  the  presence  of  an  arbitrary  distribution  of  pair  densities. 
However,  calculations  based  on  detailed  assumptions  regarding  the  pair  distribution  would  require 
a  correspondingly  detailed  microscopic  sample  characterization. 
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In  the  absence  of  a  detailed  sample  characterization,  we  can  estimate  the  magnitude  of  the 
pair  contributions  to  dielectric  enhancement  by  assuming  that  all  of  the  pairs  have  the  same 
(average)  separation  parameter.  Then  only  the  product  of  the  pair  density  and  the  incremental 
polarizability  enters  the  generalized  Clausius-Mossotti  relation  (9)  via  Equation  (8).  In  general,  a 
range  of  densities  and  incremental  polarizabilities  yield  the  same  enhancement. 

Isolated  particle  pairs  can  give  rise  to  a  substantial  enhancement  of  the  long  wavelength 
dielectric  constant  of  a  medium  of  metal  spheres.  Large  enhancements  require  large  numbers  of 
close  pain.  For  spheres  in  contact,  the  incremental  polarizability  is  of  the  same  order  of  magnitude 
as  that  of  a  single  isolated  particle.  Whatever  the  actual  pair  distribution,  only  very  close  pairs 
make  a  significant  contribution  re  dielectric  constant  enhancement. 


IV.3.7  Further  Experimental  Comparisons  with  Models 

In  the  preceding  section  we  have  seen  that  pair  interactions  alone  can  produce  permittivity 
enhancements  (above  the  Clausius-Mossotti  values)  as  large  as  those  often  observed  in  random 
metal-particle  composites.  The  example  shown  in  Figure  31  is  striking  because  of  the  ideal 
materials  of  that  experiment  [44];  i.e„  both  phases  liquid  and  perfect  sphericity  guaranteed  for  the 
metal  particles.  It  would  be  foolish  to  put  heavy  weight  on  one  or  the  other  of  the  two  limiting 
calculated  models,  inasmuch  as  the  assumption  of  Isolated  pairs  quickly  breaks  down  with 
increasing  packing.  The  apparent  trend  toward  convergence  between  experiment  and  model  is 
success  enough  in  an  otherwise  almost  intractable  problem.  For  recent  samples  of  the  difficulty  in 
treating  the  theoretical  problem  with  rigor,  see  Reference  43. 

Subsequent  to  Doyle's  report  we  encountered  another  ideal  experiment,  in  the  chemical 
engineering  literature  on  electrical  conductance  [46].  Beds  of  solid  resin  spheres  (0.2  to  1.0  mm 
diameter)  were  fluidized  by  aqueous  solutions  of  sodium  chloride.  By  varying  resin  mateiial  as 
well  as  saline  concentration,  a  wide  range  of  intrinsic  conductivity  ratios  could  be  attained. 
Volume  fractions  ranging  from  as  low  as  0.04  to  as  high  as  0.6  were  achieved.  Conductivity 
measurements  were  carried  out  at  1.6  kHz  using  specially  prepared  electrodes  in  the  fluidized  bed. 
Turner's  work  provides  a  nice  overlap  as  well  an  extension  of  the  range  studied  by  Guillien  [44], 

In  Figure  32  we  show  the  superposition  of  the  reduced  permittivity  and  reduced 
conductance  data  in  comparison  with  several  model  calculations.  First  one  should  observe  the 
rather  smooth  overlap  between  the  two  sets  of  experiments.  Together  these  data  sets  provide  an 
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Figure  32,  Reduced  permittivity  or  conductance  data  from  nveral  ideal  experiment! 

(Reference!  44  and  44)  compared  with  several  multipoie  model!  and  the  classical 
dipole  model. 
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empirical  reference  solution  to  the  problem  of  an  artificial  dielectric  composed  of  non-percolating 
(“coated”)  randomly  arranged*  spheres.  Second,  we  note  again  how  well  the  oversimplified  pair 
interaction  model  (n,  ■  plpc),  comes  close  to  describing  these  experimental  observations.  This  is 
in  strong  contrast  to  the  quasi-empirical,  quasi-theoretical  model  of  McKenzie  et  al.  [33],  which 
accounted  for  higher  order  multipoles  but  was  keyed  by  scaling  to  the  packing  limit  at  pc  ■  0.63. 
The  base  reference  curve,  the  dipolar-only  model  of  Maxwell  and  Clausius-Mossotti,  is  also 
shown.  Our  delight  with  this  modeling  success  led  us  to  describe  it  at  an  APS  meeting  [47]. 


78 


V.  MEASUREMENT  RESULTS:  PERMITTIVITY  AND  PERMEABILITY 

In  this  section  we  present  the  experimental  results  obtained  during  this  program  and 
compare  them  with  the  model  behaviors  discussed  in  Secdon  IV,  some  of  which  had  some  other 
experimental  support  and  some  of  which  had  none.  We  begin  with  a  description  of  the  several 
experimental  equipments  employed,  along  with  a  few  examples  of  typical  measurement  outputs 
(Section  V,l).  We  continue  with  presentation  of  results  for  the  polymeric  binder  artificial 
dielectrics  based  on  Nt-Cr  alloy  powders  (Section  V.2).  This  is  followed  with  some  permittivity 
results  from  polymeric  binder  composites  obtained  separately  on  an  IR&D  program  using  another 
coated  alloy  powder.  This  series  provides  an  instructive  comparison  with  the  Ni-Cr  series,  a 
comparison  which  turns  out  to  be  key  in  validation  of  the  model  behaviors.  (Section  V.3).  Finally 
we  present  our  limited  results  on  composites  fabricated  with  the  alternative  inorganic  binder 
(Section  V.4). 


V.l  Measurement  Equipment  and  Examples 

High  frequency  measuring  equipment  for  this  study  involved  a  considerable  variety  of 
approaches,  as  well  as  experiments  in  three  laboratories. 


V.  1.1  Evaluation  at  Centimeter  Wavelengths;  /  <  20  GHz 

For  the  centimeter  wavelength  range  /  <  20  GHz,  several  “Network  Analyzers”  were 
employed:  a  somewhat  antiquated  one  at  GE-CRD,  based  on  the  rather  old  (mid  1960's)  HP 
(Hewlett-Packard)  8410  system,  and  two  versions  of  the  new  HP8510  family  of  equipment. 

Measurements  of  the  constitutive  parameters,  complex  and  complex  er  were  first  made 

in  centimeter  wavelength  range  -6-18  GHz  at  GE-CRD  with  equipment  based  on  a  Hewlett- 
Packard  8410  system.  (In  principle  the  range  of  frequency  covered  can  extend  lower,  but  in  our 
practice  it  did  not  go  lower  than  6  GHz,)  The  major  building  blocks  of  the  system  were  the  8620A 
Sweep  Generator,  the  8410A  Network  Analyzer,  the  841 1 A  Harmonic  Frequency  Converter,  the 
8412A  Phase-Magnitude  Display,  the  8746B  S-Parameter  Test  Set,  and  the  90S  A  Coaxial  Sliding 
Load  for  calibration  runs. 
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Analog-to-digital  conveners,  together  with  a  data  logger  and  a  T.I.  Terminal  made  it 
possible  for  the  Network  Analyzer  to  interface  with  a  CRD  mainframe  computer.  Any  given 
frequency  band  of  the  8620A  Sweep  Generator  could  be  swept  with  a  step  size  of  0.1, 0.2, 0.5, 1, 
2,  or  5%.  Steps  of  1%  were  chosen  for  this  characterization  work.  This  enabled  many  data  points 
to  be  taken  so  that  spurious  points  were  of  less  significance  and  confidence  could  be  generated  that 
just  cannot  be  obtained  when  widely  separated  spot  frequency  measurements  are  taken.  This 
“overview”  approach  was  facilitated  by  a  screen  display  of  the  calculated  p  and  e  values  from 
which  display  a  hard  copy  could  be  obtained.  Smoothing  by  visual  inspection  and/or  a  drawing 
curve  or  straight  edge  completed  the  process.  Alternatively,  a  tabulation  of  the  measured  values 
could  be  obtained. 

Several  methods  of  measuring  the  complex  permeability  and  permittivity  are  suitable  for 
use  with  a  network  analyzer.  Of  these,  the  reflection  technique  (Sn  scattering  coefficient) 
employing  first  a  single-thickness,  then  a  double-thickness  sample  placed  in  front  of  an  RF  short 
was  the  technique  used  at  this  stage. 

Since  the  S-parameter  terminals  of  the  network  analyzer  have  precision  APC-7-type  coaxial 
connectors,  the  sample  holder  was  made  to  interface  as  smoothly  as  possible  with  this  type  of 
connector.  The  precision  outer  (0.2756-in.,  7  mm)  and  inner  (0.1 197  in.)  conductors  are  available 
commercially  and  were  used  to  construct  the  coaxial  line  sample  holder.  The  dimensions  of  a 
toroidal  sample  are  0.275  in.  O.D.,  0.120  in.  I.D.,  and  about  0.040  in.  thick.  (The  sample 
thicknesses  of  the  pair  used  for  a  particular  set  must  be  identical).  An  attempt  was  made  to  permit 
no  more  than  0.001  inch  clearance  between  the  sample  and  its  holder.  In  Figures  33  and  34  we 
show  examples  of  these  measurements  for  p'  and  fi”  in  the  6-12  GHz  band  for  sample  ONR  #2;  p 
-  0.41,  sphere  diameter  in  the  ranges  20  p.  m  to  37  p  m.  In  Figure  35,  we  show  real  permittivity 
data  for  three  loadings  of  this  same  size,  again  over  the  same  frequency  range,  p  «  0.2  to  •  0.4. 
The  figures  show  the  plotted  data  and  the  smooth  curves  used  for  subsequent  evaluation. 

During  a  visit  to  NRL  in  November  1985  to  discuss  millimeter  wavelength  experiments, 
we  brought  our  7  mm  diameter  coaxial  toroids  for  possible  confirmatory  measurements  on  their 
newer  HP  Series  8510A  analyzer.  Through  the  cooperation  of  Dr.  F.J.  Rachford  and  with  the 
help  of  S.  Browning,  it  proved  possible  to  measure  all  eight  members  of  the  series  (Table  3)  from 
zero-loading  p,  up  to  p  *  04.  In  most  cases,  more  than  one  toroid  was  measured.  The  range  of 
measurements  was  250  MHz  to  18  GHz,  but  owing  to  thickness  effects,  the  results  were  most 
useful  between  2  GHz  and  18  GHz. 
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Figure  93.  Relative  rial  permeability  measurements,  /i' ,  vs.  frequencyi  range  642  OHr, 

Sample  ONR  #2,  p  ■  0.411,  sphere  diameter  20  <  4  <  37  iim.  Measurements  on 
HP0410  system. 


FREQUENCY  GHz 


Figure  34.  Relative  loss  permeability,  vs ./»  same  equipment  and  frequency  range  as  In 

Figure  33. 
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Flgur*  35.  Rslativs  real  psrmittivlty,  t',  vs.  frtqutneyi  rings  642  OH*.  Samplas  ONR 
#2,3,4.  MtasursmtnU  on  HP  8410  systtm. 


An  example  of  the  quality  of  this  work  is  shown  in  Figure  36  for  the  sample  ONR  #2,  with 
a  volume  loading  p  ■  0.41  of  alloy  particles  in  the  size  range  20  p  m  <  d  <  37  p  m.  The  whole 
series  was  measured  in  about  four  hours!  The  overall  quantity  and  quality  of  the  results  inspired 
us  to  present  a  brief  report  at  the  March  1986  meeting  of  the  American  Physical  Society  [48]  with 
focus  on  the  permittivity  behavior  as  a  function  of  loading.  These  results  will  be  shown  in  Section 
V.2. 


The  performance  of  the  NRL  instrument  left  such  a  profound  impression  on  us  that  we 
were  able  to  convince  our  management  to  order  similar  equipment  in  mid- 1987.  (This  was  forced 
by  the  terminal  illness  and  breakdown  of  our  older  8410  system.)  The  system  obtained  is  the 
HP8510B  comprised  of  the  HP8510B  network  analyzers,  the  HP8341B  synthesized  sweeper  (10 
MHz  -  20  GHz),  the  HP8515B  S-parameter  test  set  and  an  HP  Model  9000  series  300  computer. 
A  coaxial  sample  9000  Series  holder  in  the  form  of  a  10  cm  APC-7  Beadless  Air  Line  was  used  in 
this  work.  The  measurement  procedure  is  described  in  an  HP  Product  Note  [49],  following  a 
concept  first  outlined  by  Nicholson  and  Ross  [50].  The  S-parameters  measured  provide 
information  on  the  complex  ratios  of  energy  passing  in  various  directions.  A  single  sample  is 
sufficient  when  complex  data  are  obtained  both  in  transmission  and  reflection.  With  phase  and 
amplitude  for  each,  we  have  enough  information  for  a  determination  of  all  four  constitutive 
parameters.  We  were  able  to  get  up  to  speed  on  this  machine  in  fairly  short  order  thanks  to  a 
software  package  we  purchased  from  Innovative  Measurements  Solutions  of  Atlanta,  GA,  with  the 
advice  and  assistance  of  our  colleague  R.E.  Latta  of  GEAE-Evendale,  OH. 

We  have  since  come  to  understand  some  mildly  deprecatory  remarks  by  Rachford  and 
Browning  about  the  quality  of  their  November  1985  data,  in  response  to  our  joyous  praise  thereof. 
Thus  we  now  understand  that  data  acquired  more  slowly,  with  repeat  data  averaging,  forward  and 
reverse  propagation  direction  averaging,  as  well  as  very  careful  attention  to  sample  fit  and 
thicknesses,  and  to  calibrations,  can  produce  even  better  results.  After  the  period  of  this  report,  we 
were  able  to  measure  one  of  the  original  samples  on  our  new  HP8510-B  Network  Analyzer. 
Typical  very  good  results  for  the  same  sample  #2  are  shown  in  Figure  37.  Unfortunately,  most  of 
the  original  samples  had  lost  their  mechanical  integrity  by  the  time  of  this  measurement,  so  the  full 
series  could  not  easily  be  repeated.  (The  polyurethane  binder  is  intended  for  laboratory  model 
composites  and  was  chosen  for  ease  of  fabrication  rather  than  long-term  durability.)  It  should  be 
stressed  that  the  showing  of  comparative  data  is  not  done  in  a  sense  of  competition,  but  rather  to 
illustrate  the  potential  of  the  modem  instrument. 
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Figure  36.  Example  or  constitutive  parameter  data  obtained  on  the  HP  8510A  network 
analyzer  at  NRL  for  sample  #2,  p  ■  0.41,  20  |im  <  d  <37  (im. 


Figure  37.  Example  of  same  data  as  in  Figure  36.  on  toroid  from  same  sample  material,  as 
obtained  with  greater  care  on  the  HIP  8510B  network  analyzer  at  GE*CRD. 
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With  data  of  the  attractive  quality  shown  in  Figure  37,  we  could  not  resist  upgrading  the 
results  already  obtained,  even  though  contract  time  and  funds  had  run  out.  A  limited  series  of  new 
samples  was  prepared,  and  results  from  these  are  presented  in  Section  V.2. 

V.1.2  Evaluation  at  Higher  Frequencies:  /  >  20  GHz 

Initially  it  was  planned  in  the  contract  proposal  that  all  measurements  above  the  centimeter 
wavelength  range  (f  >  18  GHz),  specifically  at  35  and  94  GHz,  would  be  carried  out  at  the  Naval 
Research  Laboratory  within  the  limits  of  their  equipment  and  scheduling.  After  our  first 
visit/discussion  at  NRL  in  February  1985,  we  were  somewhat  dismayed  by  their  requirement  of  a 
5  in.  x  5  in.  seamless  square  sample  for  the  35  GHz  measurement  in  the  arch  method.  Fortunately 
we  were  able  to  contract  with  our  colleague,  J.O.  Hanson,  at  GE  Re-Entry  Systems  Department  in 
Philadelphia  for  fixed  frequency  (35  GH2)  measurements  using  a  slotted  line  technique  [51]  with  a 
circular  waveguide  cavity  (see  Section  n.7.5).  Measurements  at  NRL  were  to  be  attempted  by 
several  techniques,  including  the  arch  method  at  rather  high  frequencies  (-90  GHz),  and 
waveguide  techniques  at  several  high-frequency  bands,  e.g.  Ka  (26-40  GHz),  Q  (40-60  GHz)  and 
W  (75-110  GHz). 

A  considerable  effort  was  expended  at  NRL  on  this  project  by  Rachford.  Satisfactory 
results  were  obtained  for  the  waveguide  samples  at  Ka  band,  as  judged  by  consistency  with  the 
lower  frequency  data  and  the  overlapping  experiments  at  35  GHz  by  Hanson.  Most  of  the 
remaining  effort  was  directed  toward  obtaining  free  space  transmission  on  an  x-ray  diffraction  table 
for  fixing  the  various  angles  needed  in  the  analysis.  The  samples  had  been  cut  to  thickness  suitable 
for  W-band.  There  were  a  few  successes  in  the  W-band  range,  but  when  these  samples  were  used 
at  Ka  or  Q  bands,  the  results  were  unsatisfactory. 

An  example  of  success  at  W-band,  the  data  for  ap  ■  0.39, 20  ^  m  <  d  <  37  fx  m  composite 
are  shown  in  Figure  38.  The  permittivity  value  of  about  20  agrees  with  Figure  35  and  36  and  the 
relative  permeability  of  about  0.6  will  turn  out  (below)  to  fit  well  with  other  results  (also  cf.  Figure 
36).  Implicit  in  the  general  fit  of  the  permittivity  is  the  condition  that  the  dielectric  behavior  of  the 
permuthane  binder  at  ~  100  GHz  is  essentially  the  same  as  the  lower  frequencies.  This  is  nicely 
confirmed  by  direct  measurements  at  NRL  on  a  polyurethane  slab  supplied  by  us.  Rachford 
reports  a  nearly  constant  value  for  e'  of  2.6  when  fitting  data  from  75  to  100  GHz.  This  is  quite 
close  to  the  lower  frequency  values  of  2.8  and  3.0  previously  obtained.  If  we  had  full  confidence 
in  the  accuracy  of  the  slightly  lower  value,  it  would  imply  an  8%  reduction  in  the  composite  £'  at 
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Figure  38. 


Permittivity  and  permeability  vs.  frequency  In  W>band  for  Nl92Crg  alloy  powder 
(20  tun  <  d  *  37  nm)  In  polyurethane,  p  *  0.39,  ai  meaiured  at  NRL. 


that  frequency  relative  to  a  low  frequency  value.  However,  while  this  spread  is  within  the  scatter 
limits  of  the  experimental  results,  such  a  decrease  is  in  line  with  the  trend  for  this  material  at  lower 
frequencies  seen  in  other  work  at  our  laboratory. 

Returning  to  the  waveguide  measurements,  several  techniques  were  employed.  The  “most 
acceptable”  results  came  from  analysis  of  network  analyzer  measurements  of  the  amplitude  and 
phase  of  the  transmitted  signal  S21  for  single  and  double  thicknesses.  When  reflection 
measurements  S11  with  and  without  a  short  circuit  backing  were  included,  there  was  considerably 
more  frequency  dependence  to  the  calculated  e  and  fx  values  than  with  the  other  experimental  set, 
as  well  as  more  than  the  “expected”  dependence,  based  on  the  lower  and  higher  frequency  data 
already  explored. 

In  Figure  39  we  show  permittivity  and  permeability  results  obtained  by  fitdng  to  the 
observed  S21  parameters  for  a  p  ■  0.42, 20jim<d<37^m  composite  sample  in  the  Ku  band. 
The  e'  value  of  24.5  is  in  excellent  agreement  with  other  results  and  the  '  behavior  will  be  seen 
as  satisfactory.  A  similar  success  was  obtained  for  the  smaller  particle,  d  <  10  /1  m,  p  »  0.36 
sample  in  this  same  frequency  band,  from  the  S21  parameters.  Deviations  from  “perfect” 
agreement  are  within  the  scatter  band  of  general  results. 


V.2  Results  -  Polymeric  Binder  Composites  with  Nl-Cr  Alloy  Powders 

In  presenting  these  experimental  results  obtained  with  Ni-Cr  powders  in  polyurethane,  we 
start  with  the  permeability  measurements.  At  this  stage  we  also  include  comparison  with  the 
modeling  of  magnetic  dipole  and  particle  size  effects  from  Section  IV.  1.  We  follow  with 
presentation  of  permittivity  measurements,  but  defer  comparison  with  models  until  Section  V.3. 

V.2. 1  High-Frequency  Permeability  Results  and  Model  Comparisons 

We  now  present  the  ensemble  of  permeability  data  superposed  on  model  calculations  of  the 
type  illustrated  in  Figures  17-20.  These  include  the  results  obtained  during  the  contract  term  plus 
the  additional  measurements  with  the  CRD-HP85 10B  system  which  we  just  couldn't  resist  doing. 
The  former  group  includes  the  frequency  and  size-dependent  results  from  the  2-18  CHz  range 
from  the  NRL  HP85 10A,  in  addition  to  our  earlier  data  obtained  from  the  outdated  HP84 10 
system,  the  just  described  Ka  band  results,  and  the  35  GHz  fixed  frequency  measurement,  and  the 
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Figure  39.  Permittivity  and  permeability  vs.  frequency  In  K.-band  for  NlpiCrg  alloy  powder 
(20  Mm  <  i  <  37  M>»)  I®  polyurethane,  p  ■  0.42,  ae  measured  by  waveguide 
techniques  at  NHL. 
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one  data  group  (p  -  0.4)  from  W  band  (75  to  1 10  GHz).  (Loss  permeability,  /t",  was  not 
measurable  in  the  NRL  experiments  above  20  GHz.)  The  gathering  together  of  all  these  results 
introduces  some  small  confusion  on  some  graphs,  since  scatter  is  inevitable.  (The  worst  is  for  p  = 
0.4  which  embraces  all  three  particle  sizes.)  Consequently  we  plot  separately  the  results  from  the 
CRD-HP8510B,  but  pair  these  graphs  with  their  counterparts  from  the  earlier  data  sets. 

We  also  confess  that  we  did  not  directly  measure  the  volume  loading,  p,  values  on  the  final 
set  as  would  have  been  obtained  using  low-temperature  saturation  magnetization  data  {>s,  at  6K), 
the  method  described  in  Section  n.7.3  herein.  Instead  for  this  set  only,  we  used  tb<-:  e'  values 
referred  to  the  ef  vs.  p  data  base  on  similar  composites  (from  the  following  subsection)  to  estimate 

p  values  and  to  verify  that  these  were  adequately  close  to  the  nominal  intended  loadings.  (The 
discussions  to  follow  on  permittivity  will  validate  this  estimation  procedure.) 

The  pair  of  graphs  of  \i\  vs.  frequency,  p  -  0.1,  data  superposed  on  model  calculations  are 
presented  as  Figures  40  and  41.  The  three  symbols  for  data  in  Figure  41  represent  three 
measurement  frequency  bands  for  which  we  usually  changed  sample  thickness  in  order  to  improve 
the  quality  of  the  results.  These  bands,  naturally  each  with  its  own  calibration  runs,  were  0.045  to 
2  GHz,  1  to  5  GHz,  and  5  to  20  GHz.  Sample  thicknesses  should,  at  most,  barely  exceed  Xm  /  4 
(see  Section  n.7.5),  at  the  highest  frequency  of  a  measurement  range.  On  the  other  hand,  a 
desirable  but  not  mandatory  minimum  thickness  would  be  AM  / 10  at  the  lowest  frequency  of  a 
range.  When  compromises  are  made  for  reasons  of  time  efficiency,  it  is  this  desired  minimum 
thickness  which  can  be  relaxed  with  only  a  modest  sacrifice  in  data  quality. 

The  graph  pain  for  p. '  vs.  frequency  continue  through  Figure  47.  Since  the  follow-up 
series  (for  the  HP8510B)  was  not  complete,  not  all  the  prior  cases  of  volume  loading  and  panicle 
diameter  (or  frequency  range)  were  duplicated.  However,  the  trends  to  agreement  as  a  function  of 
frequency  and/or  particle  size  are  qualitatively  confirmed  in  the  earlier  data  sets,  and  rather  nearly 
quantitatively  confirmed  without  spurious  inflection  points  in  the  follow-up  >  lata  sets.  In  Figure 
47,  the  upward  offset  ford  <  10  ^m,  /-range  from  5  to  20  GHz,  is  attributed  to  a  small  error  in 
the  calibration.  If  the  offset  is  ignored,  as  it  should  be  in  such  a  case,  the  predicted  curvature  is 
well  reproduced. 

The  graph  pain  for  the  loss  or  imaginary  part  of  the  permeability  follow  in  Figures  48- 
55.  The  relative  smoothness  of  the  follow-up  measurement  curves  from  the  HP8510B,  in 
comparison  to  erratic  waviness  in  the  earlier  sets  of  data,  is  a  further  testimonial  to  the  instrumental 
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Figure  42.  Comparison  of  experiment  and  model  calculations:  vs,  frequency,  p  •  0.2; 

material  system  of  Figure  40.  Measurements  at  GE  and  NRL. 


Figure  43.  Follow-up  measurements  on  HP  SS10B:  /i'  vs.  frequency,  superposed  on  model 
calculation:  composite  equivalent  to  counterpart  In  Figure  42. 


91 


Figure  44.  Comparison  of  experiment  and  model  calculations:  vs.  frequency,  p  ■  0.3; 

material  system  of  Figure  40.  Measurements  at  GE  and  NRL. 


Figure  45.  Follow-up  measurements  on  HP  8510B:  n‘r  vs.  frequency,  superposed  on  model 
calculation:  composite  equivalent  to  counterpart  in  Figure  44. 
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FREQUENCY,  GHz 

Figure  46.  Comparison  of  experiment  and  model  calculations:  /*'  vs.  frequency,  p 
material  system  of  Figure  40.  Measurements  at  OE  and  NRL, 
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Figure  47.  Follow-up  measurements  on  HP  8510B:  /i'  vs.  frequency,  superposed  on  model 
calculation)  composite  equivalent  to  counterpart  in  Figure  46. 
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Figure  48.  Comparison  of  experiment  and  model  calculations:  jur"  vs.  frequency,  p 
material  system  of  Figure  40.  Measurements  at  GE  and  NRL. 
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Figure  SO.  Comparison  of  experiment  and  model  calculations:  n?  v®*  frequency,  p  ■  0.2; 
material  system  of  Figure  40.  Measurements  at  GE  and  NRL. 
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Figure  51,  Follow«up  measurement*  on  HP  8510B:  va.  frequency,  superposed  on  model 

calculation;  composite  equivalent  to  counterpart  In  Figure  50. 


Figure  52.  Comparison  of  experiment  and  model  calculations:  vs.  frequency,  p  ■  0.3; 

material  system  of  Figure  40.  Measurements  at  GE  and  NRL. 
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Figure  53.  Follow*up  measurements  on  HP  S510B:  vi.  frequency,  superposed  on  model 

calculation;  composite  equivalent  to  counterpart  In  Figure  52. 
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Figure  54.  Companion  of  experiment  and  model  calculation!:  n"  v».  frequency,  p 
material  system  of  Figure  40.  Meaiurements  at  GE  and  NHL. 
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Figure  55.  Follow'up  measurement!  on  HP  8510B:  m,"  v».  frequency,  superposed  on  model 
calculation;  composite  equivalent  to  counterpart  In  Figure  54. 
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capability  when  extra  care  is  exercised.  In  these  graphs  we  have  “pushed”  to  extract  lower 
frequency  results  as  fx"  decreases  to  rather  small  values.  The  manufacturer's  discussion  of 
instrumentation  errors  [49]  warns  the  user  to  expect  gradual  loss  of  accuracy  as  the  loss  tangent 

or  tH  /  e')  goes  below  0.1.  In  these  loss  measurements,  most  results  fall  into  this 

cautionary  region.  For  this  reason  we  ate  reluctant  to  attach  significance  to  the  modest  inflections 
displayed  in  many  of  these  fx”  vs.  f  curves.  Indeed,  we  can  probably  attribute  the  tendency 
toward  flattening  at  a  low  plateau  in  the  range  of  ju"  equal  to  0.02  or  0.03  to  these  instrumentation 
errors  inherent  in  the  method  of  measurement.  Another  factor  of  possible  relevance  is  the  role  of 
deviations  from  particle  sphericity,  a  point  to  which  we  return  in  the  succeeding  sections. 

Overall,  with  the  ensemble  of  data  from  up  to  six  different  measurement  experiments 
covering  from  about  1  OHz  to  about  100  GHz,  we  feel  very  secure  with  the  agreement  between  the 
induced  complex  permeability  model  and  experiment. 


V.2.2  High-Frequency  Permittivity  Resuits 

Our  early  measurements  of  the  constitutive  parameters  in  the  6- 12  GHz  range  on  the 
outdated  HP8410  system  yielded  reasonable  complex  permeability  values  but  some  rather  poor 
results  for  permittivity.  The  uphill  struggle  in  this  search  is  somewhat  painfully  documented  in  the 
series  of  three  annual  technical  reports  under  this  contract  [32-34]. 

By  fortunate  contrast,  the  slotted-line,  single  frequency  measurements  at  33  GHz  produced 
results  that  generally  conformed  to  more  reasonable  behavior.  Even  here,  some  data  sets  were 
unacceptable  when  physically  meaningless  negative  values  (of  non-negligible  magnitude)  turned  up 
in  computing  complex  n  and  e.  Values  of  e*  are  expected  to  be  negligibly  small  in  this  system, 
while  values  of  /i"  may  be  small,  although  perhaps  not  very  small.  Thus  the  larger,  real  quantities 
may  be  nearly  correct,  even  when  the  imaginary  components  are  less  reliable  but  not  outlandishly 
so. 


A  word  here  about  the  dielectric  behavior  of  a  metal.  Following  appropriate  models  [24, 
23, 30],  the  frequency  region  for  a  significant  electric  loss  is  governed  by  the  intrinsic  electron 
plasma  frequency  of  the  metal  and  the  electron  relaxation  time.  These  correspond  to  frequencies  in 
the  ultraviolet  and  visible  ranges,  respectively,  clearly  far  away  from  the  microwave  and  millimeter 
frequency  ranges  of  present  interest 
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For  £* ,  repeated  measurements  at  33  OHz  provided  from  two  to  six  closely  grouped 
values  for  each  of  the  seven  composite  samples  of  the  series.  When  we  combine  these  OE-RSD 
results  with  the  NRL-HP8510A  results  at  the  selected  frequency  of  10  GHz,  we  obtain  the  curve 
of  Figure  36.  The  agreement  between  the  two  sets  is  generally  quite  good,  reinforcing  the 
conclusions  of  rather  small  frequency  dependence  and  independence  of  particle  size.  The  former 
conclusion  is  also  rather  well  supported  from  Figure  37  over  the  near-decade  of  frequencies,  2  to 
18  GHz.  In  this  same  connection,  there  is  also  a  good  accord  with  results  on  highly  loaded 
samples  obtained  at  NRL  at  the  higher  frequencies  (Ka  and  W  bands)  described  in  Section  V.  1 .2 
and  Figures  38  and  39.  It  is  also  worth  noting  in  Figure  36  that  there  is  a  smooth  departure  from 
the  binder-only  value  with  increasing  p, 

We  also  see  in  Figure  36  a  marked  upward  deviation  of  the  data  curve  from  the  Maxwell* 
Mitoff  approximation.  We  have  come  to  expect  this  from  earlier  work  (Figures  26  and  27)  and 
will  examine  it  more  closely  in  the  follov  lug  s  action, 


V.3  Results  -  Polymeric  Binder  Composites  with  Another  Alloy  Powder 

In  this  section  we  present  a  few  results  from  another  study,  an  IR&D  investigation  carried 
out  at  GE-CRD  at  the  behest  of  GE-AE  (Aircraft  Engine)  in  Evendale,  OH.  It  involves  a 
composite  system  of  coated  magnetic  alloy  particles  in  the  same  polyurethane  binder  as  employed 
for  the  Ni-Cr  alloy  powder,  and  strongly  resembles  the  B-series  discussed  in  Section  IV.2.4  herein 
(from  ref.  1).  We  call  this  set  the  C-series,  and  we  discuss  its  permittivity  behavior  in  relation  both 
to  the  improved  permittivity  modeling  (Section  IV.3)  and  the  results  obtained  for  Ni-Cr  alloy 
powder  composites  (Section  2.2).  The  comparison  is  quite  useful  for  the  goals  of  the  present 
study. 


V.  3.1  High  Frequency  Permittivity  Results  and  Comparisons 

The  alloy  powder  of  the  C-series  was  prepared  at  an  outside  vendor  by  gas-atomization 
(cf.  Section  II.4).  Subsequently  it  was  divided  by  air-classification  at  GE-AE  into  a  number  of 
close-lying  size  groups.  The  sizes  were  evaluated  by  three  particle  size  counters,  Vortec,  Coulter, 
and  Microtrac,  as  well  as  by  scanning  electron  micrographs  (SEM).  We  focus  here  on  data  from  a 
group  of  five  size  fractions  called  30R,  35R,  40R,  45R,  and  30R,  Their  average  sizes  (averaged 
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Figure  56.  Relative  permittivity)  e',vi.  volume  loading)  p.  for  NlpgCrg  alloy  .at  varloue 
particle  eliee  la  polyurethane  composites.  Data  at  10  GHx  front  NRL  and  at  35 
GHs  from  GB*RSD. 
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group  of  five  size  fractions  called  30R,  35R,  40R,  45R,  and  50R.  Their  average  sizes  (averaged 
from  the  three  counters  noted)  range  from  11.0  mm  to  14.6  mm.  Figure  57  shows  SEM  photos  of 
the  powder  fraction  designed  35R,  average  diameter  1 1.8  mm.  The  sizing  appears  to  be  rather 
narrow  and  the  particles  are  quite  spherical. 

In  most  cases,  the  various  sized  powders  were  prepared  in  composites  with  about  seven 
different  loadings.  Thus  this  group  of  five  sizes  represents  about  35  samples.  The  permittivity 
results,  at  10  GHz,  as  measured  on  the  HP8510B  at  GE-CRD,  are  shown  in  Figure  58,  as  a 
function  of  volume  loading,  p.  It  might  seem  that  the  solid  line  on  the  graph  is  drawn  to  fit  the 
ensemble  of  data.  In  actual  fact  this  predrawn  line  is  the  Guillien-Tumer  Empirical  Curve,  rescaled 
to  e*  -  2.8,  i.e.,  a  line  through  the  data  points  of  the  ideal  experiments  of  References  44  and  46, 
shown  on  a  reduced  basis  in  Figure  32.  The  agreement  is  remarkably  good,  meaning  that  we  have 
three  widely  different  experimental  families  that  follow  this  same  behavior.  Data  that  fall  below  the 
line  are  often  correlated  with  samples  showing  significant  porosity.  We  do  not  have  a  glib 
explanation  for  the  few  data  values  above  the  line,  but  can  remark  that  in  re-examinadon  of  some 
such  cases,  a  density  calculation  error,  or  a  sample  heterogeneity  was  found. 

We  conclude  that  the  Guillien-Tumer  Empirical  Curve  is  both  a  reliable  guide  forjudging 
good  quality  in  systems  which  confirm  to  its  basic  assumptions  (i.e.,  spherical,  non-contacting 
metallic  particles)  and  a  usable  guide  for  estimating  volume  loading  from  permittivity  data  in 
conforming  systems.  We  also  note  that  the  fairly  close  proximity  of  the  Guillien-Tumer  data  and 
these  C-series  data  to  the  Doyle  Pair-Interaction  model  line,  shown  in  Figure  32,  reinforces  the 
qualitative  significance  of  the  model  line. 

Now  we  can  return  to  the  Ni-Cr  alloy  comjmsites  in  order  to  compare  their  e'  vs.  p  data 
with  these  other  criteria.  In  Figure  59  we  present  the  isemble  of  these  data  (from  Figure  56)  in 
comparison  with  the  Guillien-Tumer  empirical  curve  (which  is  equivalent  to  the  C-series  data 
curve).  There  remains  a  marked  upward  deviation  of  the  Ni-Cr  data  above  about  p  ■  0.2,  although 
the  initial  slope  behavior  (p  near  zero)  conforms  with  both  the  Guillien-Tumer  curve  and  the 
MaxwelVClausius-Mossotti  model.  We  are  led  back  to  the  SEM  views  in  Figure  8  of  the  Ni-Cr 
alloy  powder  which  showed  that  non-spherical  deviations  were  fairly  common.  Their  effect  on  the 
polarizability,  as  for  the  pairs  of  Figure  29,  will  be  still  another  enhancement  factor. 

Finally,  the  apparent  reliability  of  the  e'  vs.  p  curve  for  approximating  the  loading  value, 
given  a  measured  e»,  in  the  reasonably  ideal  case,  suggests  that  one  can  extend  the  procedure  even 
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Figure  57.  SEM  photo  of  the  alloy  powder  fraction  designated  35R. 
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Figure  89.  Reduced  permittivity  e'/ejvi.  volume  loading  (e£  ■  2.8)  for  Nl-Cr  alloy  powder 

compoitteii  compared  with  the  Gullllen-Turner  Empirical  Curve  and  the 
Maxwell/Clauelus'Moiiottl  Model  Curve 


to  non-ideal  cases.  This  is  indeed  what  was  done  in  Section  V.2.1  with  the  follow-on  set  of  Ni-Cr 
composites  which  were  prepared  in  order  to  refine  the  permeability  behavior  using  the 
newest  equipment. 

From  these  data  comparisons,  despite  the  deficiencies  of  the  Ni-Cr  powder,  we  find  that 
the  agreement  of  the  permittivity  behavior  of  the  related  (but  more  nearly  ideal)  C-series  of  alloy 
powders  signals  a  satisfying  level  of  understanding  as  regards  the  permittivity  of  random  metal* 
particle,  non-percolating,  artificial  dielectrics. 


V.4  Results  -  Ni-Cr  Powder  Composites  with  Alternative  Inorganic  Binder 

The  background  for  this  part  of  the  study,  that  of  demonstrating  the  role  of  binder 
permittivity  on  composite  properties,  was  introduced  in  Section  HI.  It  was  followed  by  a 
description  of  the  search  for  a  suitable  material  and  processing  technique  to  use  with  our  Ni-Cr 
alloy  powder.  We  settled  on  a  low-melting  lead-solder  glass,  from  Owens-Illinois,  designated 
SO-7.  After  various  failed  attempts  to  machine  satisfactory  toroids  for  coaxial  line  microwave 
characterization,  we  decided  to  confine  our  tests  to  the  simpler  solid  disks  needed  for  the  slotted- 
line  technique  used  at  the  single  frequency,  35  OHz,  by  our  colleague,  J.  Hanson,  at  OE-RSD 
(Philadelphia), 

A  new  series  of  samples  was  fabricated  for  these  measurements.  Their  initial 
characterization  parameters  are  shown  in  Table  9.  Although  the  least  aggressive  fabrication 
procedure  of  the  prior  work  was  that  of  sintering  the  N192  Crs  alloy  powder  with  the  low  melting 
glass,  the  resulting  composites  did  not  have  a  high  density  and  were  still  rather  weak  for 
subsequent  machining.  Therefore,  a  mild  hot  pressing  at  relatively  low  temperature  was  used  for 
most  of  the  new  samples.  Based  on  the  early  work,  this  would  likely  produce  slight  deformation 
of  the  alloy  powder,  which  could  influence  the  electromagnetic  results. 


V.4.1  High-Frequency  Permittivity  Measurements 

Measurements  on  solid  disks  were  made  at  35  OHz  for  these  samples  and  for  the  binder 
glass.  The  latter  exhibited  a  dielectric  constant,  e' ,  of  8.2,  almost  three  times  that  of  the  organic 
polymer  polyurethane,  Thus  a  clear  binder  permittivity  effect  should  be  seen.  In  Figure  60  we 
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Table  9 


thermomagnetic  EVALUATION  PARAMETERS  OF  INORGANIC  COMPOSITES 
USED  FOR  ELECTROMAGNETIC  STUDY 


Material:  SG*7  Lead-Solder  Glass  with  RS-67  Ron  of  NlnCri,  10<d<20  tm 


Sample 

Identity 

Nominal 
Alloy  v/o 

Curie 

Temp 

K 

Sat-Mag 

*»(«K) 

emu/g 

Specific 

Grav. 

Measured 
VoL  Load 
P 

SG-7-lA(hp) 

10% 

131  ±24 

3.47 

4.50 

.115 

SG-7-1B  (hp) 

10% 

149±24 

3.22 

4.47 

.106 

SG-7-2A  (s) 

20% 

141  ±23 

5.48 

4.67 

.196 

SG-7-3B  (hp) 

30% 

151±28 

7.70 

5J28 

.298 

Notes: 

hp  ■  processing  by  hot  pressing  around  600  *C  and  300  psi 
s  «  processing  by  sintering  only  at  S3S  *C  for  20  min. 

Sample  2B  broke  during  machining. 

Sample  3A  exhibited  gross  conductivity. 
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Figure  (0.  Relative  permittivity,  e',  vi.  volume  loading,  p,  for  NICr  alloy,  10  (im  <  d  < 
20  |Am,  in  lead«solder  glau  binder,  measured  at  35  GHz. 
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show  the  results  for  permittivity  versus  volume  loading,  e'  vs.  p,  for  this  series.  It  should  be 
compared  with  Figure  56,  and  the  qualitative  si.nilarity  is  obvious. 

Next  we  consider  the  comparison  on  a  quantitative  level,  presented  in  Figure  61.  This 
graph  takes  e' t  e'k  for  each  binder  where  e'b  is  2.8  for  the  polymer  and  8.2  for  the  glass.  At  low 
loadings,  up  to  20%  by  volume,  the  scaled  agreement  is  excellent.  There  is  appreciable  deviation 
at  higher  p-values, 

A  possible  reason  for  the  deviation  is  the  enhanced  polarizability  of  non-spherical  particles 
such  as  might  have  occurred  during  the  hot-pressing  preparation.  It  should  also  be  recalled  that  the 
powder  preparation  run  RS-67  already  yielded  a  modest  share  of  non-spherical  "pickle-like" 
particles,  as  was  shown  in  Figure  8.  A  comparative  scaled  plot  of  polymer  composites  of  this 
alloy  with  the  Guillien-Tumer  Empirical  curve  (which  matches  the  quite  spherical  C-series  alloy 
powder  composite  permittivity)  was  shown  in  Figure  59.  The  Ni-Cr  alloy  series  (from  RS-67) 
permitdvity  data  had  a  considerable  upward  deviation  from  the  behavior  of  the  C-series,  even 
without  any  possible  deformation  of  the  particles  in  composite  preparation.  In  hindsight,  we  now 
attribute  that  deviation  between  the  two  experimental  systems  to  the  non-spherical  nature  of  some 
of  the  Ni-Cr  alloy  RS-67  particles.  (We  note  in  passing  that  Figure  61  also  includes  35  GHz  data 
from  the  K„-band  sweeps  at  NRL  on  the  two  polymeric  composites  described  earlier.) 

To  sum  up  this  section,  the  data  of  Figure  61  do  confirm  the  role  of  binder  permittivity  as  a 
scaling  parameter  in  the  dielectric  behavior  of  artificial  dielectric  composites.  We  remarked  earlier 
that  this  has  an  impact  on  application  design.  It  also  could  be  relevant  for  quality  control  in 
composites.  If  there  were  significant  porosity  in  a  particular  composite  preparation,  it  could  be 
treated  according  to  mixing  rules  as  a  diluent  to  the  binder  [21 , 52],  which  in  turn  could  lower  the 
final  composite  permittivity. 


V.4.2  High-Frequency  Permeability  Measurements 


These  inorganic  composite  samples  measured  at  35  GHz  also  displayed  relevant 
permeability  results,  along  the  lines  previously  modelled  for  the  polymeric  binder  composites. 
Inasmuch  as  both  binders  are  non-conducting  and  non-magnetic  (/i*  *  1,^'=  0),  there  are  no 
changes  introduced  in  the  expected  permeability  behavior  by  changing  binder.  Since  a  single 
particle  size  group  was  used  for  this  series,  10  /u  m  <  d  <  20  p  m,  we  make  our  comparison 
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Figure  61.  Scaled  companion  of  permittivity  behavior,  at  35  GHi,  for  artificial 

dielectrics  with  two  different  binder!.  Nl-Cr  particles  in  polyurethane  ■  2.8] 
and  In  lcad>solder  glass  [e[  »  8.2]. 


between  experiment  and  calculation  with  the  expected  behavior  for  d  »  20  jt  m.  The  comparison  is 
shown  in  Figure  62  for  both  real  and  imaginary  parts  of  the  permeability.  The  bars  on  the 
experimental  represent  the  full  spread  of  several  data  points  obtained  from  various  disks  of  the 
same  nominal  (or  identical)  loadings.  The  agreement  is  fully  satisfactory,  a  tribute  to  the  quality  of 
the  measurement  method  and  to  the  adequacy  of  the  model. 
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Figure  <2.  Comparlaon  of  experiment  and  model  eaiculationsi  and  va.  frequency, 

varloua  loading*)  for  NiCr  alloy  powder,  10  <  4  <  20  nm,  in  lead«aolder 

glaaa.  Data  at  3 5  OH*  only,  from  QE»RSD. 


VI.  CONCLUSIONS 


This  research  has  been  concerned  with  the  development  of  an  interesting  class  of  artificial 
dielectrics  and  a  deeper  understanding  of  its  high-frequency  electromagnetic  properties.  This  Final 
Technical  Report  -  Part  A  deals  with  essentially  non-magnetic  composite  dielectrics.  Part  B, 
which  contains  classified  material,  treats  magneto-dielectric  composites. 

The  key  ingredient  of  the  artificial  dielectric  was  a  powder  of  spherical  Ni-Cr  alloy  particles 
of  ultimate  composition  Ni92  Crg  prepared  by  gas  atomization  with  liquid  substrate  quenching. 
The  latter  step  of  the  preparation  produced  an  insulating  oxide  coating  on  each  particle.  The 
composites  fabricated  with  this  metal  powder  and  polyurethane  binder  (or  with  an  alternative 
inorganic,  glassy  binder)  were  ferromagnedc  only  below  room  temperature  Tc  ■  160K  This 
feature  allowed  for  easy  determination  of  volume  loading,  using  low  temperature  saturation 
magnetization  measurements.  Before  mixing  with  the  binder,  the  alloy  particles  were  classified  by 
size  at  <10  mm,  +10-20  n  m,  and  +20-37  n  m. 

While  the  original  focus  was  intended  to  be  on  permitdvity  behavior,  the  realization  that 
induced  magnedc  dipole  (eddy  current)  effects  would  be  equally  important  led  to  the  decision  to 
pay  close  attention  to  powder  size.  Calculations  of  complex  magnedc  permeability,  adapted  from 
the  literature,  showed  that  the  dominant  parameter  is  the  ratio  of  particle  radius  to  skin  depth.  The 
sizes  chosen  and  frequencies  of  study  (centimeter  to  millimeter  wavelength  ranges)  covered  the 
most  interesting  behavior  anticipated. 

An  exploration  of  theoretical  developments  and  history  in  the  dielectric  behavior  of 
heterogeneous  media  successfully  addressed  one  of  our  original  questions,  i.e.,  why  our  prior 
(and  present)  material  systems  appeared  to  follow  the  simple  Maxwell  (Maxwell-Gamett,  Clausius 
and  Mossotti)  theoretical  model.  The  inherently  isolated  metallic  particles  represent  a  correlated 
system  even  when  loaded  randomly  in  a  binder,  in  contrast  to  the  percolation-prone  Bruggeman 
model  exemplified  in  many  experimental  metal-insulator  composites.  Modem  theories  which 
incorporate  higher  order  electric  multipole  effects  and  geometrical  packing  limits  have  been 
validated  previously  only  (to  our  knowledge)  by  ordered  lattice-like  arrays  of  spherical  particles. 
Our  random  composites  exhibited  similar  properties,  albeit  with  deviations  probably  reflecting  on 
their  microstructure. 

The  present  study  was  successful  in  improving  the  modeling  of  permittivity  behavior  of 
randomly  packed  artificial  dielectrics.  This  was  a  twin  success,  i.e.,  on  the  empirical  and  on  the 
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"theoretical”  fronts.  The  empirical  success  rests  on  the  discovery  in  the  literature  of  several  rather 
different  sets  of  results  from  ideal  experiments  involving  truly  spherical  particles,  well  randomized 
and  non-percolating.  Some  of  our  own  data  sets  strongly  resemble  this  ideal  behavior.  On  the 
theoretical  side,  we  developed,  through  the  work  of  our  consultant,  W.T.  Doyle,  a  rather  simple 
pair-interaction  model  which  is  a  naive  oversimplification.  However,  it  shows  the  importance  of 
local  clustering  and  higher-order  multipole  effects.  This  pair  interaction  description  comes  closer 
to  describing  the  experimental  results  than  any  other  model  known  to  us.  This  convergence,  even 
if  not  taken  seriously,  gives  one  a  feeling  of  confidence  and  predictive  control  in  considering  the 
dielectric  behavior. 

Measurements  of  the  electromagnetic  constitutive  parameters  (permittivity  and  permeability) 
for  Ni-Cr  alloy  powder/binder  composites  were  carried  out  in  three  laboratories.  Centimeter 
wavelength  (2-20  GHz)  results  were  obtained  at  GE-CRD  and  at  the  Naval  Research  Laboratory, 
Higher  frequency  experiments  were  carried  out  at  GE  Re-Entry  Systems  Department  (Philadelphia) 
at  35  GHz,  and  at  NRL  for  Ka-band  (26-40  GHz)  and  W-band  (75-110  GHz).  Also  included  in 
our  experimental  analysis  were  data  obtained  in  another  study  under  internal  sponsorship  using  a 
different  alloy  powder  (C-series  data).  Its  particles  have  a  higher  quality  morphology;  i.e.,  there 
were  fewer  deviations  from  sphericity  than  for  the  Ni-Cr  alloy. 

From  the  hierarchy  of  steadily  improving  experimental  results,  the  induced  permeability 
data  are  in  good  quantitative  agreement  with  the  calculated  model  over  the  parameter  ranges 
examined,  i.e.,  diameters  from  about  5  /im  to  35  n  m  and  frequencies  from  about  1  GHz  to  100 
GHz. 


The  permitdvity  data  are  Independent  of  particle  size,  and  change  but  little  over  the 
frequency  range  covered.  The  latter  small  effect  is  probably  due  to  properties  of  the  binder.  The 
main  interests  in  the  permittivity  focused  primarily  on  its  dependence  on  volume  loading,  and 
secondarily  on  the  scaling  properties  with  change  of  binder.  The  C-series  permittivity  data  are  in 
remarkably  good  agreement  with  the  empirical  curve  obtained  from  ideal  experiments  in  the 
literature,  our  so-called  Guillien-Tumer  Empirical  Curve.  The  experimental  curve  has  values  only 
slightly  higher  than  those  predicted  by  the  oversimplified  pair-interaction  model  of  Doyle.  As  for 
the  Ni-Cr  powder  composite  peimittivity,  a  residual  upward  deviation  of  e'  vs.  p  is  attributed  to 
non-spherical  departures  in  particle  morphology.  These  are  visible  in  scanning  electron 
micrographs, 
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From  the  experiments  with  an  alternative  inorganic,  glassy  binder,  we  confirm  the  scaling 
normalization  law  that  e*  /  £*  (c  ■  composite,  b  «  binder)  puts  different  systems  on  an  equivalent 
footing.  Indeed  this  was  built-in  to  early  theory  and  had  been  confirmed  many  times  before. 

This  research  has  demonstrated  that  the  microwave  electromagnetic  behavior  of  non¬ 
percolating,  artificial  dielectric  composites  is  well  understood.  Random  arrays  have  their  own 
unique  properties,  as  distinct  from  ordered  arrays.  Whereas  several  reports  in  the  late  1970's  and 
early  1980's  [24,25]  had  pointed  to  giant  deviations  from  electromagnetic  predictions  for  far- 
infrared  absorption  of  submicron  metal  particles,  more  recent  studies  [18]  showed  qualitatively  that 
inadequate  inter-particle  dispersion  was  the  problem.  Our  work  puts  this  explanation  on  a 
quantitative  basis,  albeit  at  modestly  lower  frequency  and  larger  particle  size.  The  present  study 
also  provides  a  firm  foundation  for  examining  the  more  difficult  metallic  magneto-dielectric 
composites  discussed  hi  Part  B  of  this  Final  Report. 
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